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FOREWORD 


This  report  was  prepared  by  Lexington  Laboratories,  Inc., 
Cambridge,  Massachusetts  under  USAF  Contract  AF33 (657) -11280 . 
This  contract  was  initiated  under  Project  No.  7360,  "The 
Chemistry  and  Physics  of  Materials",  Task  No.  736001,  "Thermo¬ 
dynamics  and  Heat  Transfer" .  The  work  was  administered  under 
the  direction  of  the  Air  Force  Materials  Laboratory,  Research 
and  Technology  Division,  Wright-Patterson  Air  Force  Base,  Ohio, 
with  Mr.  R.  J.  Prezecki  acting  as  project  engineer. 

Robert  C.  Folweiler  is  Project  Engineer  on  the  contract 
for  Lexington  Laboratories,  Inc.  Contributors  to  the  project 
include  W.  B.  Fraser,  D.  A.  Hill,  W.  D.  Kingery,  W.  J.  Mallio, 
and  D.  W.  Moore. 


ABSTRACT 


Thermal  radiation  from  refractory  oxides  is  a  volume, 
rather  than  a  surface  phenomenon  in  the  wavelength  range 
0.5-6  microns.  Emissivity  values  have  been  calculated  for 
specially  prepared  AI2O2,  MgO,  Si02  and  SrTiO^  samples  from 
appropriate  reflectivities,  absorption  coefficients  and 
scattering  coefficients  derived  from  material  character¬ 
istics  (index  of  refraction,  single  crystal  transmissivity, 
pore  size  and  concentration) . 

Experimental  measurements  of  emissivity  for  these 
samples  have  been  made  at  tempei ^ture  up  to  1200°C  and 
at  wavelengths  from  one  to  fifteen  microns.  Calculated 
and  measured  values  are  in  good  agreement  for  ^^2^3 
SiOj.  Absorption  coefficients  of  the  MgO  and  SrTiO^  samples 
were  greater  than  the  single  crystal  values,  leading  to 
higher  measured  (than  calculated)  emissivities.  Samples 
of  commercial  alumina  had  higher  emissivities  than  the 
specially  prepared  high  purity  samples. 

The  Hamaker  equations  for  non- isothermal  emissivity, 
which  linear i.- a  the  temperature  gradient,  have  been  com¬ 
pared  with  the  exact  solution  by  numerical  integration. 
Results  are  in  good  c.greement  for  all  cares  of  practical 
interest  for  application  to  ceramic  oxides.  Solutions  to 
the  Heunaker  equations  have  been  derived  for  new  boundary 
conditions. 

Evaluation  of  the  rotating  blackbody  slot  system  in¬ 
dicates  it  gives  results  in  good  agreement  with  an  external 
blackbody.  An  apparatus  for  measuring  the  effective  emis¬ 
sivity  of  a  non- isothermal  sample  l\as  been  constructed  and 
tested . 
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This  technical  dociimentary  report  has  been  reviewed  and 
is  approved. 


Assistant  Chief 
Materials  Physics  Division 
AF  Materials  Laboratory 
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Net  radiative  flux  defined  by  wF  =  I~J 

Special  function  defined  in  Ref.  5,  Appendix  I 

Absorptive  coefficient  used  by  Ryde 

Special  function  defined  in  Ref.  5,  Appendix  I 

Net  heat  flux 

Radiant  energy  flux  in  the  direction  of  the  positive 
x-axis 

Radiant  energy  flux  in  the  direction  of  the  positive 
x-axis  at  X  =  0 

Specific  intensity  of  radiation 

Specific  intensity  of  forward  radiation  at  x  =  0 

Parallel  forward  flux 

Parallel  flux  incident  on  surface  of  slab 
Average  specific  energy 

Specific  intensity  of  refracted  radiation 
Specific  intensity  of  radiation  scattered  once 
Specific  intensity  of  radiation  not  scattered 
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Emission  coefficient  for  once  scattered  radiation 
Radiant  energy  flxix  in  the  direction  of  the  negative 
x-axis 

Radiant  energy  flux  in  the  direction  of  the  negative 
x>axis  at  X  -  0 

k  Lattice  thermal  conductivity 

K  Scattering  factor  (non-dimensional) 

Jl  Path  length  in  direction  0 

m  Reciprocal  index  of  refraction;  m  =  1/n 

n  Index  of  refraction 

M  Number  of  pores  per  unit  volume 

p  Dimensionless  parameter  p  =  4Trr  im-l| /km 

P  Vol\ime  fraction  porosity 

r  Pore  radius 


s  Backward  scattering  coefficient  for  diffuse  radiation 

(scattering  factor  in  Hamaker  equations) 
s'  Backward  scattering  coefficient  for  parallel  radiation 

S  Scattering  coefficient  for  parallel  radiation 


S 


diffuse 


Scattering  coefficient  for  diffuse  radiation 


t  Dimensionless  parameter  t  =  (o'  +  S)x 

t'  Dimensionless  parameter  t'  =  (»  +  S)x' 

T  Absolute  temperature  - 

U  Material  constant  defined  by  Eq.  2.69 

V  Material  constant  defined  by  Eq.  2.70 

X  Distance  along  x-axis 

O'  Absorption  coefficient 

3  Material  constant;  3  =  o/(a  +  2s) 

3  Material  constant;  3  =  c  /(a  +  2s) 

o  o  o 


y  Euler's  constant;  y  =  0.'j772156 
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Constant  denominator  defined  in  Eq.  4.23 
Emissivity 
Normal  eiaissivity 
Hemispherical  emissivity 

Angle  measured  from  x-cocis  within  the  sample 

Angle  iraasured  from  x>cucis  outside  the  sample 

Material  constant  defined  by  Hamedcer 

_  2^ 

”  ka 


Wavelength  of  radiation  in  free  space 
Reflectivity 

Reflectivity  of  diffuse  energy  emerging  from  the 
sample 

Reflectivity  of  diffuse  energy  incident  on  the  sample 
Reflectivity  for  normal  incidence 


Material  constant  for  non- isothermal  case 


Material  constant 


=  '>ra(a+2s) 


Stefan-Boltzmcun  constant 


Azimuthal  angle 
Spectral  radiant  energy 
Solid  angle 
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1.0  INTRCDUCTICN 


The  evaluation  of  ceramic  oxides  for  high* temperature 
applications  and  the  dasign  of  high* temperature  structures 
require  a  precise  knowledge  of  the  thermal  gradients  in  and 
radiant  heat  transfer  to  and  from  solid  bodies.  Energy  trans¬ 
fer  through  the^e  systems  is  by  two  modes:  radiation  and 
lattice  conduction.  Interactions  between  these  heat  trans¬ 
fer  modes  can  cause  large  deviations  frcxn  linear  temperature 
gradients  near  surfaces. 

Emissivity  has  been  widely  used  as  a  material  constant 
characterizing  what  is  often  considered  as  a  surface  phenomena. 
However,  for  oxide  materials  which  are  partially  transparent 
in  the  important  0.5-8  micron  wavelength  range,  radiant  energy 
emission  is  a  volume  process.  As  a  result  there  is  no  simple 
"emissivity"  for  a  material  like  alumina,  but  the  radiant 
energy  emitted  depends  on  temperature  gradients  and  also  on 
material  characteristics  such  as  absorption  and  scattering 
coefficients.  Thus,  the  n<^menclature  often  used  for  opaque 
solids  with  "emittance"  taken  as  a  sample  characteristic, 
while  "emissivity"  is  a  material  constant  is  not  applicable. 

We  v;ili  use  the  term  "emissivity"  for  the  ratio  of  energy 
emitted  from  an  isothermal  sample  to  the  energy  emitted  from 
an  ideal  black  body  at  the  same  temperature,  fully  recognizing 
that  this  is  a  sample  characteristic  rather  than  a  material 
constant. 

In  the  present  investigation  we  have  been  mainly  con¬ 
cerned  with  evaluating  the  material  characteristics  for  well 
defined  systems,  developing  a  method  for  calculating  the 
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isothermal  emissivity  from  these  characteristics,  and  compar¬ 
ing  the  calculated  results  with  direct  experimental  measure¬ 
ments.  This  is  an  ambitious  approach,  but  if  successful,  it 
allows  us  (a)  to  estimate  the  thermal  emissivity  from  material 
characterization  rather  than  having  to  make  expensive  and  dif¬ 
ficult  mt-^sureraents  on  every  sample,  and  (b)  to  evaluate  the 
effect  of  changing  material  characteristics  on  thermal  emis¬ 
sivity  behavior  and  thus  provide  a  rational  basis  for  new 
material  development. 

1  2 

In  addition  we  have  evaluated  the  basic  Ryde  -  Hamaker 
equations  for  calculating  and  applying  isothermal  and  non- 
isothcrmal  emissivities ,  extended  the  Hamaker  relation  to 
new  boundary  conditions,  and  developed  a  technique  and  appara¬ 
tus  for  carrying  out  controlled  measurements  of  non-isothermal 
emissivity. 
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2.0  CALCULATION  OF  ISOTHERMAL  EMISSIVITy  PROM  MATERIAL 
CHARACTERIZATT  ON 

The  material  factors  which  enter  the  expression  for 
isothermal  emissivity  are  the  reflectivities  for  entering 
and  leaving  radiation  and  the  material  constant, 

3  =  >/  a/(a  +  2s) 

o 

where  a  and  s  are  the  absorption  and  scattering  coefficients 
discussed  subsequently.  The  emissivity  for  a  thick  sample 
is  given  by; 


(1  -  e)  = 


(1-p^)^  +  2B^a-p^^)  +  8^^(i+p.)^ 


where  is  the  reflectance  of  diffuse  radiant  energy  emergent 
from  the  solid  and  p  is  the  reflectance  of  diffuse  radiant 
energy  incident  on  the  solid.  Thus  calculation  of  emissivity 
requires  evaluation  of  p^,  p^,  and  3^. 

As  discussed  by  Folweiler^  the  ratio  a/s  can  be  measured 
independently  of  the  emissivity  by  measuring  diffuse  trans¬ 
missivity  of  thin  samples.  Unfortunately,  the  measurements 
are  not  very  precise.  For  example,  two  measurements  of  trans¬ 
missivity  for  beryllia  gave  values  within  +3%.  However,  calcu¬ 
lations  of  =  n/  a  (a  +  2s)  gave  values  different  by  a  factor 
of  two  and  the  a/s  ratios  varied  by  a  factor  of  five;  this 
leads  to  a  difference  in  calculated  emissivity  of  a  little 
less  than  a  factor  of  two.  The  magnitude  of  these  errors 
results  from  the  fact  that  the  calculations  depend  on  a  rela- 
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tively  small  difference  between  large  measured  values,  and  is 
inherent  in  the  form  of  the  mathematical  relationships  between 
variables. 

2 . 1  Reflectivity 

The  total  reflectivity  for  diffuse  radiation  incident 

on  a  surface  depends  on  the  index  of  refraction  of  the  solid 

1  4 

as  discussed  by  Ryde  and  Ryde  and  Cooper.  Values  for  a  range 
of  refractive  indices  are  listed  in  Table  2.1. 

Table  2.1 

Total  Reflectivity  for  Diffuse  Radiation  and  for  Radiation  at 

4 

Normal  Incidence  (After  Ryde  and  Cooper  ) . 

Total  Reflectivity 

Index  of  Refraction  normal  incidence  diffuse  radiation 


1.00 

0 

0 

1.1 

0.0023 

0.026 

1.15 

0.0049 

0.035 

1.2 

0.0083 

0.045 

1.25 

0.012 

0.053 

1.3 

0.017 

0.061 

1.35 

0.022 

0.069 

1.4 

0.028 

0.077 

1.45 

0.033 

0.085 

1.5 

0.040 

0.092 

1.55 

0.047 

0.100 

1.6 

0.053 

0.107 

1.65 

0.060 

C.114 

1.7 

0.067 

0.121 

1.3 

0.082 

0.134 

1.9 

0.096 

0.146 

4 


From  the  values  of  p  given  in  2,1,  we  can  obtain 
the  reflectivity  of  emergent  radiation,  p^,  from  the  relation: 


1  - 


2 

n 


(2.2) 


2.2  The  Absorption  Coefficient,  a 

The  absorption  coefficient,  a,  in  the  Hamaker  rela¬ 
tion  refers  to  the  absorption  of  energy  from  diffuse  radiation 
as  it  passes  through  a  body  containing  scattering  centers.  In 
contrast,  the  intrinsic  absorption  coefficient,  Of ,  is  measiured 
by  and  refers  to  the  absorption  of  energy  from  parallel  radia¬ 
tion  in  a  body  without  scattering  centers. 

Thus,  application  of  measured  values  of  a  to  calcula¬ 
tion  of  must  take  into  account  the  larger  path  length  for 
diffuse  radiation  as  it  passes  through  a  thin  layer  of  thick¬ 
ness  dx.  Let  us  consider  the  specific  intensity  of  radiation, 
I' ,  where  9  measures  the  angle  frcxa  the  x  direction  (see 
Pig.  2.1).  We  consider  situations  in  which  I'  varies  with 
angle  9  from  the  x  direction,  but  is  independent  of  the  azi¬ 
muthal  angle  (i.e.,  symmetrical  about  the  x  axis);  that  is, 

I'  =  I'(x,9).  Consider  first  an  intensity  field  in  free  space 
specified  by 


I’(x,9)  =  I'q  for  0  ^  9  <  Tr/2,  (2.3) 

and  I'  (x,9)  =0  for  ir/2  <  Q  <  ir . 

This  field  is  isotropic  over  the  left  and  right 
hemispheres  separately,  and  the  flux  associated  with  this 
field  is 

TT  F  =  I  I'cos  9  dcD  (2.4) 

ID 

where  the  integration  is  over  all  solid  angles. 
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Then, 

2tt 

TT 

irp  « 

II  • — » 
o 

I 

9*0 

1*  (x,9)  cos  9  sin  9  d9  d<^  * 

tt/2 

2ir 

I 

0*0 

I' 

o 

cos  9  sin  9  d9  *  ¥l '  . 

o 

(2.5) 

So,  ?  *  I'  . 

o 

Now  suppose  that  instead  of  free  space  we  specify 
this  radiation  field  for  0<9<ir/2  on  the  plcme  x  =  0  in  a 
material  with  absorption  coefficient  ^  per  unit  length.  At 
X  -  d  we  specify  that  no  radiation  is  incident  on  the  material 
from  the  right.  Then,  the  intensity  I'  must  satisfy 

^  =  -  or  I'  (x,e)  (2,6) 

where  i  measures  path  length  in  the  direction  8.  clearly, 
i  =  X  sec  8.  So, 


4~  sec  0  1'  , 

CDC 

I '  (x,0)  s=  I '  (0,9)  exp  (-tt  X  sec  9]  0  <  6  <  ^  , 

cuid  I'(x,9)  =  I'(d,9)exp  I-o(x-d)sec  9]  ^9<ir,  (2.7) 

But  we  specify  I'  (0,9)  *  I’^for  0  <  9  <  Tr/2, 

and  I'  (d,0)  =0  for  ■w/2  <  9  <  ¥. 
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I'(x,0)  =  I'q  exp  [-  »  X  sec  0]for  O<0<tt/2 

and  l'(x,0)  =0  for  n/2<0^TT. 

The  flux  is  given  by 

7rp(x)  =  I  l'(x,0)  cos  0  do). 

n/2 

'irF(x)  =  27r  J  I'  exp(-ttx  sec0)  cos0  sin  0  d0 
o  o 

—  1  rj  sef*  6 

Since  cos  0  =  (sec  0)  ;  - _  =  sin  0  d0,  and 

(sec  0) 


F(x)  =21 


» 

« 


exp{-ot  X  sec  0)  -  =  2  I'  I 

1  (sec  0)  ° 


where  E*  is  the  exponential  integral  of  order  3  given  by 

^  5 

Chandrasekhar,  Appendix  I,  (p.  373). 

For  a  small  argument,  has  the  power  series 

expansion 

£^(2)  =  1/2  -  z  +  1/2  (-7- In  z  +  3/2)  ^  +  .  . 

where  y  is  Euler's  constant  (see  Van  de  Hulst,  p.  244). 
P(x)  =  I'q  (1-2Q'x  +  0  (a^x^))  . 


,  (2.8) 
(2.9) 

(2.10) 

.  (2.11) 

S3  (orx) 
(2.12) 

.  .  (2. 

Thus, 

(2.14) 


So  for  optically  thin  slabs ,  we  have 


(2.15) 


F(x)  -  F(0) 
P(0) 


2a  X. 


That  is,  the  equivalent  absorption  coefficient  for  diffuse 
radiation  in  a  thin  slab  is 

a  *  2a.  (2.16) 


When  there  are  scattering  centers  present  the  analysis 
is  complicated  by  the  fact  that,  as  a  result  of  scattering,  some 
path  lengths  are  increased  and  some  path  lengths  are  decreased. 
Thus  it  is  not  obvious  whether  scattering  will  increase  or  de¬ 
crease  a.  If  we  consider  a  slab  of  material  between  x  =  0  and 
X  =  d,  let  the  material  contain  a  random  distribution  of  scat¬ 
tering  centers  that  scatter  isotropically  with  a  scattering 
coefficient  S  per  unit  path  length.  We  consider  the  situation 
in  which  the  intensity  just  inside  the  left  boundary  (at  x  *  0+) 
is  given  by 

I'(0+,  0)  *  I*  for  0  <  0  <  Tr/2,  (2.17) 

and  just  inside  the  right  botmdary  (at  x  =  d-)  is  given  by 

I ' (d- ,  0)  =0  for  7r/2  <  0  <  tt  , 

Let  I '  be  the  intensity  of  radiation  which  has 
ns 

not  been  scattered  at  all,  and  I*  the  intensity  of  once 

s 

scattered  radiation.  The  only  source  of  I '  is  at  the  bound- 
ary  x  =  0.  In  the  material  is  decreased  by  two  processes. 

Some  is  absorbed  by  the  material;  this  process  is  characterized 
by  the  absorption  coefficient,  a.  The  other  process  that 
decreases  is  scattering  (this  energy  goes  into  I'^)  which 
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Pig.  2. 


Diagram  of  geometry  of  radiation  in  a  thin 
slice  of  absorbing  material. 
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NJ  IW 


may  be  characterized  by  the  scattering  coefficient  S.  Taking 
both  absorption  and  scattering  into  account, 


dl' 

d(x  sec  6) 


{a  +  s)I' 


ns 


(2.18) 


and  the  solution  satisfying  the  boundary  conditions  is 

I'  =1'  exp  (-(O'  +  S)  X  sec0]  forO  <9  <  ir/l, 
ns  o  “  “ 

and  I'  (x,9)  *  0  for  tt/2  <  0  <  n.  (2.19) 


Now  the  amount  of  radiation  scattered  in  all  direc¬ 
tions  per  unit  volxome  of  material  from  I'  (x,0)  is  SI'  (x,0)  . 

ns  ns 

Scattering  is  isotropic,  so  the  amount  scattered  per  unit 

volume  per  unit  solid  angle  is  (S/4Tr)l'  (x,0)  . 

ns 

Integrating  this  over  all  directions  of  incoming 

radiation  gives  the  emission  coefficient  for  once  scattered 

radiation,  is  defined  as  the  amount  of  radiant  energy 

emitted  per  unit  volume  per  unit  time  per  unit  solid  angle 

(due  to  scattering  of  I '  ) .  So  we  have 

^  ns 

n  2ti 


I 


I' 


ns 


(x,e) 


I 

0=0 


I  ’ 


ns 


(x,0) 


sin  0  d0  d0. 

(2.20) 


ti/2 

r 

J 

0 


n 


.1 


=  1  J 


I '  (x,9)  sin  9  d0 

ns 


I'  exp  [  -  (Q'+S)x  sec  0]  sin  0  d0. 
o 


(2.21) 
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Let  M.  -  cos  @. 


I 


exp  [  -(«+S)  ”]  d4 
M' 


[(or+s)x] 

$ 


(2.22) 


where  E.  is  the  exponential  integral  of  order  2. 

^  5 

l*he  equation  of  transfer  (see  Chandrasekhar,  p.  6) 
for  once  scattered  radiation  is 


dl '  SI ' 

sec  e^dx  =  -  +  -2^=2  (2.23) 


This  may  be  written  as 


dl 


(«+S)  sec  0  dx 


=  I 


s  2(Q'+S)  o  2 


E-  [(fl'+S)x] 


(2.24) 


For  0  <  0  <  7r/2  the  boundary  condition  is  I '  (O+,0)  =0,  and 

s 

for  7r/2  <  0  <  TT  the  boundary  condition  is  I'  (d-,0)  *  0.  The 

s 

solution  to  the  equation  which  satisfies  this  boundary  condition 
is 


I 


SI  ' 

I'  (x,0)  =  ^  J  E  [  (Q'+S)x']  exp  [-(o+S)  sec0(x-x')  ]sec0  dx' 

s  ^  ^  z 

for  O<0<^  .  (2.25) 


Let  t'  =  (Q'+S)x'  and  m-  =  cos  0.  Then, 


SI '  ^ 

l'^(x,0)  =  [-4>'+S)sec  9  x] 


SI '  sec  0 

2  -  ®XP  [“(«+S)sec’0  x]  F2(t,ii) 


J  E2  (*■-')  exp  (t'/ii)dt' 


0 


(2.26) 


(2,27) 


for  0  <  0  <  Tr/2. 
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where  t  »  (0'‘i’S)x  and  is  defined  in  Appendix  1  of  Chandrasekhar.^ 
He  gives 


p2(t,^L)  +  exp  (t/ji)  E^Ct)  -  1],  (2.28) 


So, 


SI' 

I'g(x,e)  ~2(<»+S)  (-t/pL)  Pj^(t,kL)+  E2(t)  -  exp  (-c/|i)]  (2.29) 

for  O  <  0  <  n/2 . 

Now  let  I'  *  I'  +1*  and  neglect  all  higher  order  scattering. 

ns  s 

Ihis  is  a  good  approximation  if  the  slab  of  material  is  optically 
thin.  We  wish  to  calculate  the  flux  of  energy  across  any  plane 
X  =  const.  This  is  given  by 


n 


'rrp(x)  = 


i  I'  (x,< 


I  I*  (X,l 


,0)  cos  0  dO)  =  2  n  J  I*  (x,0)  cos  0  sin  0  d0  (2.30) 
60  0 
TT 


Trp(x)  =  2tt  J  [I'  (x,0)  +1'  (x,6)l 

^  ns  s 


cos  0  sin  0  d0 


(2.31) 


Now  we  know  I'  (x,0)  =1'  exp  [-(ff+S)x  sec©]  for  0  <  0  <  Tr/2 
AS  o 

(2.32) 


and  I'  (x,0)  =0  for  Tr/2  <  0  <  tt. 
ns  "■ 


(2.33) 


Also  I'  (x,0)  is  known  for  0  <  0  <  Tr/2.  But  we  also  need  to 
s  ”  ~ 

know  i’O'r  Tr/2  <  0  ^  tt  to  do  the  complete  flux  integral. 

We  have  not  solved  for  I'^  for  thes*.  angles  and  we  do  not  need 

to  if  we  restrict  ourselves  to  calculating  the  flxix  at  x  =  d, 

where  the  boundary  condition  is  I'  (d,0)  =0,  Tr/2  <  6  <  w, 

s 
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That  is,  no  scattered  radiation  comes  from  the  right.  Using 


this  fact,  we  have 

T7/2 


irF(d)  *2Tr 


J  II'  (d,9)  +  I'  (d,0)]  cos  6  sin  6  d0  (2.34) 
ns  s 


TTSI '  P 

.  2irl'  E,t(«+S)di  J  “P  l-(“+S)d/n)  t(«+S)d,M.)d4. 

°  ^  ^0  (2.35) 

Thus  the  integral  is  recognizable  as  the  function  ^2,2^^^^^^^^ 

given  by  Chandrasekhar^  in  Appendix  1. 
so, 

irp(d)  =  2ir  !•  [E3  ( («+Sld)  +  ^ ,„®s)  ^  ([<»+S]d)]  (2.36) 


Van  de  Hulst  gives  power  series  expansion  ^2,2  ^ 

Appendix  1  of  his  1948  paper  in  Astrophys.  J. 


xn 


E2(x)  =  1/2  -  X  +  1/2  (-7  -In  X  +  3/2)  x  +  .  .  .  (2.37) 


G',  (x)  =  x+  (7+lnx  -  3/2)  x  + 

2,  f  2 


(2.38) 


where  x  >  0,  and  y  is  Euler's  constant. 

So,  for  (ff+S)d  «  1  we  have. 


TrF(d)  =  2Trl^' [*1  -  (ff^-S)d  +  0((ft+S)^d^)  + 


[  (or+S)d  +  0[(a+S)^d^]  )] 


2(a+S) 

*  2  2 
■n-F(d)  =  2Trl  '  -  (or+s/2)d  +  0  (  {a+S)  d  }] 


o  '•2 


(2.39) 

(2.40) 
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-  1  -  (20  +  s)d  +  0  l(o+s)  Vl 


(2.41) 


So, 


This  result  has  em  easy  physical  interpretation. 

For  an  optically  thin  slab  of  thickness  d,  a  fraction  2  (ar+s)d 
of  the  incident  flux  is  cd>sorbed  or  scattered.  Of  the  fraction 
2dS  which  is  scattered,  half  goes  into  backward  flux  and  half 
into  forward  flux,  since  scattering  is  isotropic.  Neglecting 
further  scattering  and  absorption,  we  see  that  a  fraction  Sd 
of  the  original  flux  is  added  to  the  forward  flux  by  scattering. 
A  fraction  2  (o'+S)d  is  removed  by  absorption  scattering  and  a 
fraction  Sd  is  added  by  the  scattering.  So  the  net  fraction 
removed  is  (2Q'+S)d«  That  is,  the  absorption  coefficient 
appropriate  to  the  Hamaker  conditions  is  not  affected  by 
isotropic  scattering  centers,  and  Eg.  2.16  is  the  appropriate 
one  to  use  for  calculating  a. 


2.3  The  Scattering  Coefficient,  s 


The  scattering  coefficient,  s,  as  used  in  the 
Hamaker  equations  refers  to  back  scattering  of  diffuse  radia¬ 
tion.  In  contrast,  the  commonly  measured  and  calculated 
scattering  coefficient,  S,  refers  to  total  scattering  for 
parallel  radiation.  Assuming,  as  Hamaker  did,  diffuse  radia¬ 
tion,  by  the  same  analysis  as  given  for  the  absorption  co¬ 
efficient  above,  ~  2S;  taking  only  the  back  scatter¬ 

ing  (assuming  isotropic  scattering)  as  shown  in  Eq.  2.41, 


s 


1 

2 


S 


diffuse 


S 


(2.42) 
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The  scattering  coefficient,  S,  is  related  to  the 
number  of  pores  per  unit  volume,  N,  and  the  pore  radius,  r, 
by  the  scattering  factor,  K,  in  the  relation 

S  =  K  N  (irr^)  (2.43) 


If  the  volume  fraction  porosity  is  P, 


(2.44) 


The  scattering  factor  K  varies  between  zero  and 
four  and  depends  on  the  pore  radius,  r,  wavelength  of  radia¬ 
tion,  'h,  and  the  relative  index  of  the  scattering  center 

7 

and  medixam  m.  Van  de  Hulst  has  summarized  solutions  for  a 
number  of  systems  in  terms  of  the  dimensionless  i<arameter, 
p  =  4Trr  jm-ll  /Am. 

8 

R.  H.  Boll,  R.  D.  Gtimprecht,  and  C.  M.  Sliepcevich 
have  calculated  the  scattering  factor  K  from  the  Hie  equa¬ 
tions  for  values  of  the  relative  refractive  index  of  0.93, 
0.9,  and  0.8  corresponding  to  our  situation  of  spherical 
pores  in  a  solid  media.  Their  results  show  that  the  excur¬ 
sions  of  K  from  an  average  value  near  2  are  less  marked  as 

9 

m  decreases.  D.  W.  Lee  and  W.  D.  Kingery  extrapolated  Boll 
et  al's  results  to  lower  values  of  m  and  compared  calculated 
and  directly  measured  scattering  coefficients  for  alumina 
and  Vycor  having  pore  sizes  in  the  range  0.7-7  microns  with 
good  results  (agreement  vdthin  about  +10%) . 

Now  a  basic  question  is  "How  isotropic  is  the 
scattering  by  pores  in  an  oxide  matrix?"  That  is,  can  we 
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use  the  available  calculations  for  the  total  scattering  co¬ 
efficient  discussed  above  to  evaluate  the  back  scattering 
appropriate  to  the  Hamaker  relation? 

The  angular  intensity  distribution  of  the  scattered 
radiation  depends  on  the  size  of  the  scattering  center,  the 
wavelength  of  radiation,  and  the  relative  refractive  index. 

For  very  small  particles  (relative  to  the  wavelength)  scatter¬ 
ing  is  nearly  isotropic.  For  very  large  particles  having  an 
index  of  refraction  greater  than  the  media  (m  >  1) ,  the  scat¬ 
tering  is  almost  completely  forward.  However,  for  the  reverse 
case,  pores  of  lower  index  embedded  in  a  media  of  higher  index, 
this  is  not  the  case.  As  shown  in  Fig.  2.2  when  m  =  0,  the 

7 

back  scactering  is  just  half  the  total  scattering.  J’or  Al^O^ 
and  MgO,  the  back  scattering  is  about  25%  for  the  total  scat¬ 
tering  . 

For  wavelengths  of  0.5-8  microns  and  pore  sizes  from 
0.5-5  microns  (experimentally  observed)  we  cover  the  gamut 
from  "very  small"  to  "very  large"  pores.  Complete  solutions 
of  the  Mie  equations  would  be  required  for  the  intermediate 
range.  Since  the  relative  indsx  enters  each  term  of  the  re¬ 
quired  summations,  there  is  no  easy  way  to  transpose  available 
solutions  to  the  case  of  scattering  pores  in  a  solid  matrix. 

Fortunately,  the  relative  values  of  forward  and 
backward  scattering  for  pores  in  a  solid  matrix  cover  a  much 
narrower  range  than  is  found  for  the  more  usual  case  of  solid 
particles  in  a  matrix  of  lower  refractive  index.  For  pores  in 
alumina,  for  example,  we  can  write: 

|K£<S<|Kf  (AljOj)  (2.45) 

That  is,  we  can  calculate  the  appropriate  back  scattering  co¬ 
efficient  to  better  precision  than  it  can  be  evaluated  from 
measurements  of  diffuse  transmissivity. 
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DmATRIX 


FRACTfON  OF  BACK  SCATTERING 


Fig.  2.2  Relative  Fraction  of  Back  Scattering  for  Pores  of 
Large  Size  in  a  Solid  Matrix. 
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It  should  be  recognized  that  this  calculation  is 
appropriate  only  for  pores  in  a  solid  matrix  —  which  is  the 
system  we  are  concerned  with.  It  is  not  applicable  to  more 
complex  multiphase  compositions  in  which  scattering  centers 
have  a  higher  index  than  the  matrix  material. 

2.4  The  Material  Constant,  S 

o _ 

The  constant,  S  =  si  a/(a+2s)  can  be  calculated 

o 

directly  from  values  of  a  and  s  determined  as  described  in 

the  previous  section.  Since  this  calculation  involves  taking 

the  square  root  of  s,  the  precision  of  calculating  8^  is 

better  than  the  evaluation  of  s.  That  is,  the  maximum  expected 

two- fold  error  in  s  leads  to  only  a  ^^2  error  in  B  . 

o 

2.5  Calculation  of  Emissivity  in  the  Wavelength  Range  0-5-8 
Microns 

In  this  range  of  wavelengths  emissivity  is  a  volume 
process,  and  calculations  have  been  made  from  Eq.  (2.1) . 

Values  of  the  reflectivities  for  each  wavelength  were  taken 
from  Table  1.  Values  of  the  appropriate  absorption  coefficient 
were  calculated  for  each  wavelength  from  measured  values  of 
the  intrinsic  absorption  coefficient  with  Eq.  (2.16) .  values 
of  the  appropriate  back  scattering  coefficient  were  calculated 
for  each  wavelength  and  pore  size  from  Eq.  (2.44)  using  values 
of  the  scattering  factor  K  taken  from  Fig.  (2.3) . 

Results  are  compared  with  directly  measured  isothermal 
emissivity  in  Section  8.0. 

2.6  Calculation  of  Emissivity  in  the  Wavelength  Range  8-15 
Microns 

At  wavelengths  above  5-8  microns,  the  absorption 
coefficient  increases  to  such  an  extent  that  emissivity  becomes 


18 


Fig.  2.3  Values  of  the  Scattering  Factor  K  for  Various  Values  of 
the  Relative  Index,  m.  When  ra  <  1. 


61 


SCATTERING  FACTOR  (  K ) 

ro 


T 


t» 
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a  sxirface  process  depending  on  the  reflectivity.  Eg.  (2.1) 
reduces  to: 

1  -  e  *  (2.46) 

Bmissivity  values  can  then  be  calculated  directly  from  the 
data  in  Table  2.1 

At  variable  wavelengths,  which  are  above  10-12  microns 
for  most  oxides  the  restrahl  frequency  is  approached,  the 
reflectivity  increases  and  the  emissivity  drops  off.  We  have 
not  analyzed  behavior  in  this  wavelength  region  in  detail. 

2.7  Normal  vs.  Hemispherical  Bmissivity 

In  the  calculations  discussed  thus  far,  we  have 
been  concerned  with  the  spectral  hemispherical  emissivity, 
while  our  experimental  measurements  have  been  for  spectral 
normal  emissivity. 

For  the  wavelength  range  8-15  microns  where  emis¬ 
sivity  is  essentially  a  surface  process,  and  Eg.  (2.46)  applies, 
the  difference  is  directly  related  to  the  difference  in  normal 
and  diffuse  reflectivity  as  given  in  Table  2.1.  Comparable 
expected  emissivity  values  are  given  in  Table  2.2. 

For  the  wavelength  range  below  about  8  microns  where 
emissivity  is  a  volume  process,  we  are  concerned  with  the  an¬ 
gular  distribution  of  the  emergent  radiation.  Radiation  of 
intensity  1'^  at  angle  9  will  be  partly  reflected  and  partly 
refracted  at  the  surface.  Assuming  isotropic  diffuse  radia¬ 
tion  just  below  the  surface  we  can  calculate  the  angular  in¬ 
tensity  distribution  of  the  emergent  radiation. 
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Table  2.2 


Expected  Normal  euid  Hemispherical  Emissivity  for  Smooth 

Surfaces  Having  Given  Values  of  the  Refractive  index  (for 

Sample  with  a  »  s)  . 

Refractive  Normal  Hemispherical 

Index  Emissivity  Emissivity 

1.1 

.998 

.974 

1.2 

.992 

.955 

1.3 

.983 

.939 

1.4 

.972 

.923 

1.5 

.960 

.908 

1.55 

.957 

.900 

1.6 

.947 

.893 

1.65 

.940 

.886 

1.7 

.933 

.879 

1.8 

.918 

.866 

1.9 

.904 

.854 

Let  dl  be  the  rate  of  transport  of  radiant  energy 
per  unit  area  in  directions  confined  to  an  element  of  solid 
angle  doo  (see  Fig.  2.4) .  dl  is  then  expressed  in  terms  of 
intensity  I ' ; 


dl  =  I'  cos  e  (2.47) 

Consider  isotropic  radiation  of  intensity  I'  incident  on 
a  plane  interface.  Of  the  energy  dl  incident  on  a  unit  area 
from  the  solid  angle 

do)  =  sin  0  d0  d0  (2.48) 
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r^dl  is  reflected  and  (l-r.)  dl  is  refracted  into  the  solid 

tj  y 

angle 


cU'=  sin  6'  d0'  d0' 


(2.49) 


where 


0  =  0« 

sin  0'=  n  sin  0 


(2.50) 


The  intensity  l'0i  of  the  refracted  radiation  is  defined  by 


(l-r^)  dl  =  I‘q,cos  0'  dttt' 


(2.51) 


which,  with  the  use  of  equations  (2.47) ,  (2.48)  and  (2.49) 
becomes 


(1-r^)  I'  cos  0  sin  0  d0  d0  =  1*^,  cos  0*  sin  0'  d0'd0 


(2.52) 


From  equation  (2.50) 


cos  9'  sin  0'  d9'  =  n  cos  0  sin  0  d0 


(2.53) 


and  using  (2.53)  in  (2.52) 


I'e,  =  (-yii' 


(2.54) 


The  reflectivity  is  given  by  the  well-known  Fresnel  equation 


_  1  r  sin^(9-0')  tan^(0-0') , 

^0  2  2  "^2 

sin  (0+9* )  tan  (0+0') 


(2.55) 
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Then  for  any  value  of  n,  we  can  calculate  the 

angular  distribution  of  the  relative  intensity  of  the  reflected 

and  transmitted  energy.  This  is  shown  for  AI2O2  (n  =  1.7)  in 

Fig.  2.5.  The  relative  values  of  the  normal  and  hemispherical 

emissivity  are  determined  by  comparing  the  normal  intensity 

tr  he  average  value  integrated  over  all  wavelengths. 

The  integrated  energy  flux  rate  per  unit  area  is 

rr/2 

I  =  J  I*  ,  cos  0*  doj*  -•  ¥  i*  3in20'  d0'  (2.56) 

y  n 

hemisphere  o 

wheie  use  has  been  made  of  equation  (2 .49)and  the  integration 
over  the  azimuthal  angle  0  has  been  carried  out  to  arrive  at 
the  right-hand  expression  in  equation  (2.56).  If  1'^,  ±s  an 
isotropic  intensity,  I' ,  the  integration  over  0'  can  be 
carried  out  to  give 


I  *  ¥  I  (2.57) 

Therefore  an  "equivalent"  isotropic  intensity  will  be  defined 
by 


tt/2 


si'1  20'  dG' 


tt/2 

=  J 


n 


(l-r^(e') ]  sin  20'  dG' 


(2.58) 


Prom  equation  (2.54) 


I 


0'=O 


I ' 


2 

n 


(1  -  p  ) 


(2.59) 
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Fig.  2.5  Distribution  of  Intensity  of  Reflected  and  Trans¬ 
mitted  Energy  of  Isotropic  Radiation  Incident  on 
an  Interface. 
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where  p  =  r^,  ^  is  the  reflectivity  for  normal  incidence 
n  w  -u 

tabulated  in  Table  2.1. 

Hence  dividing  (2.59)  by  (2.58) 


I ' 

9'=Q 


tt/2 

J (l-r^)  sin  29'  dG' 
o 


(2.60) 


The  j.^i,cegral  in  the  denominator  may  be  evaluated  by 
changing  the  variable  of  integration  from  0'  to  0  according 
to  equation  (2.50)  and  noting  that  the  range  of  integration 
will  then  be  from  0  =  0  to  9  =  arc. sin  (1/n) .  Thus; 


tt/2 


r 

o 


(l-rg) 


sin  29'd0'=  n^ 


arc . sin ( 1/n) 

(1-r^)  sin  20  d0 

0 


arc. sin (1/n) 

==  1  -  n^  J  ^9 

o 


(2.61) 


This  integral  has  been  evaluated  by  Walsh, and  the 

3 

result  1.S  also  recorded  by  Polweiler  (p  18)  .  The  result  is 
finally 


it/2 

J (1-rg) sin  20 'd0' 


1  (n~l)  (3n+l)  n^  (n^-1)  ^ 

^  6  (n  Kl)  ^  (n^+1;  ^ 


log 


inzH 

(n+1) 


2n^ (n^+2n-l)  -  8n^ (n^+1) _ 

2  4  2  4  2 

n  +1)  (n  -1)  (n  +1)  (n  -1) 


=  1  - 


^o 


(2.62) 
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p  is  the  reflectivity  for  diffuse  radiation  tabulated  in 
Table  2.1. 

Hence  equation  (2.60)  may  be  written 


(2.63) 


That  is,  for  isotropic  radiation  emergent  frtm  a 
surface,  the  ratio  between  the  normal  and  hemispherical  emis 
sivities  is  given  by  the  relative  value  of  (1  -  p^)  and 
(1  -  p^) .  Values  for  different  refractive  indices  cure  given 
in  Table  2.3. 
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Table  2.3 

Ratio  of  the  Normal  to  Hemispherical  Emissivities  for  Varying 
Index  of  Refraction. 


n 

1-p 

^  n 

1-p 

^o 

^o  h 

1.0 

1.000 

1.000 

1.000 

1.2 

.992 

.955 

1.038 

1.3 

.983 

.939 

1.047 

1.4 

.972 

.923 

1.053 

1.5 

.960 

.908 

1.057 

1.6 

.947 

.893 

1.060 

1.7 

.933 

.879 

1.061 

1.8 

.918 

.866 

1.060 

1.9 

.904 

.854 

1.058 
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2.8  Reflectivity  for  Parallel  Radiant  Flux  at  Normal 
Incidence 

If  a  beam  of  parallel  radiation  is  incident  normally 
on  a  slab,  it  is  necessary  to  extend  the  calculation  for ^dif¬ 
fuse  incident  light  in  Section  3.1  of  Polweiler's  report. 

The  necessary  calculations  have  been  done  by  Ryde.^  Some 
of  Ryde's  equations  are  reproduced  here  to  have  them  available 
in  the  notation  of  this  report.  It  is  necessary  to  define 
three  new  scattering  coefficients,  f,  f ' ,  and  s' ,  where  these 
are  a  forward  scattering  coefficient  for  the  diffuse  radiation 
and  forward  and  backward  scattering  coefficients  for  the 
parallel  radiation  respectively. 

The  equations  for  the  forward  and  backward  diffuse 
fl\ix  are  then: 


j  T 

=  f  I  -  (a  +  s)  I  +  sJ 
dx  P 


(2.64) 


dJ 

T—  =  s'  I  -  (a  +  s)  J  +  3l 
dx  P 


(2.65) 


and  for  the  parallel  flux  I 


dl 


,  =  (a  +  s'  +  f )  I  =  g'  I 

dx  p  ^  p 


(2.66) 


The  solutions  to  these  equatxons  obtained  by  Ryde 


are 


I  =  e 
P 


-g  'x 


(2.67) 


ax  -ax 

I  =  A(1-8  )e  °  +  B(1+8  )  e  ° 

o  o 


-  Ue 


-g'x 


(2.68) 
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ax  -ax  , 

J  =  A (1+0  )e  °  +  B(1-8  )  e  °  -  Ve”^  ^ 

o  o 


(2.69) 


whore 


:2af '  +  (s+f )  (s'+f 
2af  +  (s'+f)^ 


(2.70) 


(s-s')  (f +8') 
[2af  +  (s'+f)^l 


(2.71) 


and  A,  B  are  to  be  determined  from  the  appropriate  boundary 
conditions,  which  are. 


I  =  (1  -  p  )I  and  I  =  p.  J  at  X  =  0; 
p  n  p.  o  L 


and  J  =  p.I^  ot  X  =  D; 
D  D 


I  is  the  parallel  flux  incident  on  the  surface  of  the  slab. 

Also,  the  assumption  is  made  that  I  Ai/  0  at  x  d  or  more 

P 

explicitly  that  g»a. 

From  the  second  and  third  boundary  conditions  we  find 


-(U-V) {(1-p^)  -  (1+p^)  3^3  exr(-a^d) 
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and  using  the  first  boundary  condition  the  reflectivity  p  of 
the  surface  is 


P  -  P„  +  r  =  Pn  *  ‘^-Pn’  <^-Pi'  ■^< 


(2.74) 


Substituting  for  J 


P  =  Pn  + 


f(U-V) {[(l-p.)  -  (l+p.)8^ 

_ Hi _ Li _ c 

([(1-p^)^  +  (l+p^)^8^^]  si 


!  y 

']  sinh  c  d  +  2  p.0  cosh  cj  d}  | 


inh  a  d  +2(l-p.  )8  cosh  a  d 
o  o  o 


(2.75) 


If  the  incident  radiation  is  diffuse,  the  reflectivity  p  may 
be  obtained  from  the  above  by  setting  f  *  f  and  s'  =  s,  in 
which  case  U  =  1  and  V  =  0;  and  by  replacing  p^  by  p^  which 
gives 


p  =  p^  +  (l-p^)  (1-p^) 


2 


(1-p^) - (1+p^) 0^  sinh  a^d  +  2p^0^  cosh  a  d 
[(1-p^)^  +  sinh  a^d  +  2(l-p^^)0^  cosh  a^d 


(2.76) 


which  is  the  same  «£)  F*^.  (20)  of  Folweiler's  report. 
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3.0  THE  HAMAKER  APPROXIMATION  FOR  CONDUCTIVE  HEAT  PLOW 

In  the  develofment  of  the  Heuncdcer  equations,  an  approxi¬ 
mation  is  introduced  for  the  conductive  part  of  the  combined 
heat  transfer  by  conduction  and  radiation.  Measurements  of 

the  effective  thermal  conductivity  of  polycrystalline  alumina 

o 

at  temperatures  near  2000  c  indicate  that  the  conductive  part 

9 

is  still  a  major  contributor,  <>ven  at  this  high  temperature. 
This  results  from  the  relacively  high  conductivity  and  large 
scattering  coefficients  in  polycrystalline  oxides. 

Since  the  utility  of  the  Hamaker  equation  depends  on  the 
validity  of  the  \inder lying  assumptions,  we  felt  it  essential 
to  critically  evaluate  the  nature  of  these  limitations. 

The  exact  equations  for  describing  the  transmission  of 
heat  by  combined  conductive  and  radiative  transfer  in  a  homog¬ 
eneous  medium  are  (for  the  one-dimensional  case  and  the  steady 
state) ; 


k  ^  -  2  a  n^  a'  t"^  +  a  (I  +  J)  =0  (3.1) 

dx 

^  =  a  n^  a'  T^  -  (a  +  s)  I  +  sJ  (3.2) 

ox 

=  (a  +  s)  J  -  sl  -  a  n^  a'  t"^  (3.3) 

uX 

where  I  and  J  are  the  radiant  fluxes  in  the  +x  and  -x  direction, 
T  is  the  absolute  temperature,  a'  is  the  Boltzmann  constant, 
n  is  the  refractive  index,  k  is  the  thermal  conductivity,  a 
is  the  absorption  coefficient  for  diffuse  radiation,  and  s 
is  the  back-scattering  coefficient  for  diffuse  radiation. 
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3.1  Direct  Numerical  Integration 

It  is  convenient  to  reduce  the  set  of  Eg.  (3.1) , 
(3.2),  and  (3.3),  to  a  single  second-order  equation.  Eg.  (3.2) 
and  Eg.  (3.3)  are  added  and  the  sum  integrated  once  to  give: 


I  +  J  =  -  (a  +  2s)  (H  X  +  kT)  +  A 


(3.4) 


in  which  A  is  a  constant  of  integration  and  H  is  the  constant 
heat  flux  defined  by: 


dT 

H^I-J-kT^ 

dx 


(3.5) 


Siibstituting  in  Eg.  3.1  we  obtain: 

k  -  2  a  n^  a'  -  a  (a  +  2s)  (Hx  +  kT)  +  a  A  =  0 
dx  (3.6) 


For  simplicity  we  will  consider  the  one-dimensional 

flow  of  heat  by  conduction  and  radiation  inside  a  medium  as 

shown  in  Fig.  3.1.  Once  the  material  constants,  k,  a,  s,  and 
2 

n  are  specified,  there  remain  four  independent  initial  con¬ 
ditions  to  be  fixed;  i.e.,  T  ,  (dT/dx)  ,  I  ,  and  J  .  For 

o  o  o  o 

various  values  of  these  initial  conditions  we  have  constructed 
a  family  of  solutions  by  numerical  integration  of  Eg,  (3.6)  . 

The  calculation  consists  of  (1)  selection  of  material 
constants  and  initial  conditions;  (2)  evaluation  of  A  and  H 
from  Eqs.  (3.4)  and  (3.5);  (3)  calculation  of  T(x)  by  numerical 
integration  of  Eg.  (3.6) ;  (4)  calculation  of  k  (dT/dx) ;  and 
(5)  calculation  of  I (x)  and  J(x)  from  the  relations: 


1  dT 

(a  +  2s)  (Hx  +  kT)  +  a  +  H  +  k  (—)  ] 
i  dx 


(3.7) 
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(3.8) 


(a  +  2s)  (Hx  +  kT)  +  C 


-  H-k(g)] 


3.2  The  Hamaker  Approximation 

Hamaker ' s  approximation  consists  of  replacing  the 

first  term  in  Eq.  (3.1)  k  (d‘^T/dx^)  Ly  d^(T^)/dx^  in 

which  is  some  average  temperature  in  the  material.  This 

replacement  makes  the  set  of  Eqs.  (3.1),  (3.2)and  (3.3;  linear 

4 

and  homogeneous  in  the  quantities  I,  J,  and  T  so  that  a  solu¬ 
tion  in  closed  form  can  be  obtained  readily. 

3.3  Numerical  Evaluations  of  the  Hamaker  Solutions 

The  Hamaker  solutions  for  the  radiant  fluxes  are 
given  in  terms  of  the  Ixamped  material  constants: 


b  =  4n  o'  T 

o 

(3.9) 

a  -  ^  +  a  (a+2s)  =  ^/  (1  +  /c) 

(3.10) 

0  =  a/(a  +  2s) 

(3.11) 

ic  =  2b/k  (a  +  2s)  =  2b0/k(r 

(3.12) 

Then,  the  radiant  energy  fluxes  and  temperature  are  given  by: 
I  =  A(1-3)  e^^  +  B(l+0)  e"^^  +  C  (ctx-B)  +  D  (3.13) 

J  =  A(l+0)e^^  +  B(l-0)e  +  C(dx+|3)+D  (3.14) 

(JX  -0X 

E  =  -A/ce  -  B/ce  +  Cax  +  D  (3.15) 

T  =  Tq  +  ^  (E_e^)  (3,16) 
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/ 


where 

J  a  (1+0)  -  E  (20a)  -  (^  )  20 
o  o  dx 

A  =  _ 2 _ 

40  a  (/c  +  1) 

dp 

-JO  (l-$)  -  (2So)  +  (;t^)  2(3 

o  o  ox 

B  = _ 2 _ 

4  8  a  (/c  +  1) 

dP 

J  2  /c  a  +  (^)  20 
o  dx 

C  =  - 2 - 

4  0a  (k  +  1) 

J  2  0  /c  a  +  E  (40  a) 
o _ o _ 

°  "  4  0a  (/c  +  1) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


The  calculation  consists  of  (1)  selection  of  material 
constants  and  initial  conditions;  (2)  evaluation  of  b,  a,  0 
and  K  from  Eq.  (3.9),  (3.10),  (3.11),  and  (3.12);  (3)  evalua¬ 
tion  of  A,  B,  C,  and  D  from  Eqs.  (3.17) ,  (3.18) ,  (3.19) ,  and 
(3.20);  and  (4)  calculation  of  I,  J,  E,  and  T  for  various  values 
of  X  from  Eqs.  (3.13) ,  (3.14) ,  (3.15) ,  and  (3.16) . 

3.4  Comparison  of  the  Exact  and  Hamakrr  Solutions 

A  comparj...on  has  been  made  between  nxmierical  inte¬ 
gration  of  the  exact  equation  as  described  in  Section  3.1  and 
numerical  evaluation  of  the  Hamaker  solution  as  described  in 
Section  3.3  for  the  system  shown  in  Fig.  3.1  for  a  range  of 
initial  conditions.  Computations  were  done  rn  an  IBM  7090 
computer  at  the  Aeronautical  Systems  Division,  Wright- 
Patterson  Air  Force  Base.  For  one  set  of  conditions  a  hand 
calculation  was  done  in  excellent  agreement  with  the  computer 
calculation  as  shown  in  Fig.  3.2. 
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Temperature 


Fig.  3.2  computer  Results  (Solid  Lines)  and  Hand 
Calculated  Results  (Points)  for  the 
Exact  Solution  of  Radiant  Energy  Heat 
Transfer . 


Heat  flux  (cal  cm’^sec" 


TEMPERATURE  CK) 


2000 


10 


0  EXACT  EQUATION 
•  HAMAKER  SOLUTION 


500 


Fiq.  3.3  Typical  results  between  exact  and  Hamaker  equations 

for  small  values  of  a  and  s.  I  *  0;  J  *  J; 

^  o 

0.1;  s  0.1; 


,dT 


k  =  0.01;  =  1,000;  (^)  *  10;  a 

T  o  dx 

n  =  7j  . 
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HEAT  FLUX  (CAL  SEC’*  CM" 


•1». 


200 


0  EXACT  SOLUTION 
•  HAMAKER  SOLUTION 


X  (MM) 

Pig.  3.4  Typical  Results  of  Comparison  between  Exact  and 
Haraaker  Equations  in  Large  Values  of  a ,  s  and 
dT 

=  0;  J  =  50;  k  =  .01;  T  =  2000; 
dx  o  o  o 

o 

dT 

(^)  =  1000,  a  =  1,  s  =  100,  n2  = 


3. 


HEAT  FI 


Table  3.1  Input  Constants  and  Comparison  of  Results  of  Exact  and  Hamaker  Equations 

for  20  Conditions 


40 


000  1.000  1.000 


Table  3.1  (Continued) 
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For  all  calculations  was  set  equal  to  zero;  physi¬ 
cally  this  corresponds  to  radiation  from  the  surfa-^e  into  a 
semi-infinite  media  with  the  same  index  as  the  test  section 
(no  surface  reflection)  and  having  perfect  transmissivity  (no 

scattering  or  absorption) .  Values  of  J  ranged  from  0.26-125 

-2  -1  ^  -1 
cal  cm  sec  ,  values  of  (dT/dx) ^  ranged  from  lO-lOOO^K  cm 

and  values  of  T  ranged  from  100-3000®K.  Material  constants 

o  -1-1  -1 

were  varied  ini  the  range:  k,  .005-05  cal  cm  sec  ®K  ; 

-1  -1  2 

a,  0. 1-1.0  cm  ;  s,  0.1-100  cm  ;  n  =3.  Some  typical  re¬ 
sults  are  shown  in  Fig.  3.3  and  Fig.  3.4.  Tabulated  input 
data  and  a  comparison  of  results  between  the  Exact  and 
Hamaker  solutions  are  given  in  Table  3.1. 

The  depth  of  the  sample  contributing  to  the  emis- 
sivity  characteristics  is  that  part  that  is  optically  "near" 
the  surface.  Thus,  even  though  calculations  were  carried  out 
to  a  greater  sample  depth  in  most  cases,  we  have  compared 
results  at  X  =  (a  +  2s)  ^  and  x  =  2(a  +  2s)  ^  in  Table  3.1. 

As  shown  there  for  the  entire  range  of  conditions  where  a 
and  s  are  greater  than  unity,  agreement  between  the  exact 
and  the  Hamaker  equations  was  within  0.1%.  For  the  calcula¬ 
tion  with  a  =  0.1  and  s  =  0.1,  agreement  was  within  about  1% 
at  X  =  (a  +  2s)  Considering  that  values  of  thermal  con¬ 

ductivity  up  to  the  highest  observed  for  oxides  (0.05  for  BeO) 
and  values  of  temperature  gradients  up  to  1000° /cm  were 
evaluated,  the  results  clearly  show  that  the  Hamaker  equation 
is  a  satisfactory  approximation  for  any  thermal  conditions 
likely  to  be  encountered  for  heat  transfer  problems  related 
to  oxide  materials. 
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4.0  EXTENSION  OF  THE  HAMAKER  EQUATIONS 


The  Hamaker  equations  for  radiant  energy  fluxes  and 
energy  distribution  are: 


I  =  A(l-3)e^^  +  B(l+0)e"‘^^  +  C(ax-j3)  +  D  (4.1) 


J  =  A(l+0)e^^  +  B(l-3)e"^^  +  C{ax+3)  +  D  (4.2) 


(JX  -(JX 

E  =  -A/ce  -  B/ce  +  Cax  +  D 


where  a  = 


■t-/Tab  + 

J  k 


a(a+2s)  =  OaJ  (l+/c) 


(4.3) 

(4.4) 


and  0  =  a/(a+2s) 


and  ic  =  2b/k(a+2s)  -  ZbB/'ka 


(4.5) 

(4.6) 


We  have  recalculated  the  values  of  the  constants  A, 

dE 

B,  C  and  D  for  boundary  conditions  when  I  ,  J  ,  (t~)  and  E 

o  o  dx  ® 

are  given  at  x  =  0  without  the  limitation  that  I  =0.  We 

o 

find; 


-20  (^) 
(dx) 


a  (1-0)  I  +  a  (1+0',  j  -  2  0aE 
o  o  c 


4  0  (7  (I  J-a;) 


(4.7) 


20  (dE)  +  CT  (1  +  0)  I  -  a  (1-0)  J  -  20OE 
3  =  (iV/0  ' 
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c 


(4.9) 


40  (^) 
(dx) 


-  2x01  +  2xaJ 

o  o 


40a  (1+x) 


D  = 


2a0x(i_+j  )  +  40aE 
_  o  o _ o 

400  (1+x) 


(4.10) 


We  have  also  calculated  values  for  the  boundary  condi¬ 
tions  I  =  0atx  =  0,T  (or  E  )  at  x  =  0 ,  T ,  or  E ,  at 
o  o  o  da 

X  =  d,  and  a  constant  for  the  heat  flux  H,  where 


,  dE 

H  =  I  -  j  -i-s 
o  o  b  dx 


(4.11) 


Inserting  the  boundary  conditions  in  equations  4.1  and 

4.2, 


=  A  (1-6)  +  B  (1+0)  -  C0  +  D 
J  =  A  (1+0)  +  B  (1-0)  +  C0  +  D 

O 

From  4.3 

dE  ,  ^  „  -crx  ,  ^ 

-r~  -  -Axoe  +  Bxoe  +  Co 
dx 

H  =  I  -  J  -'r(-  A.co  +  BXO  +  CO) 
o  o  b 

From  4.12,  4.13,  and  4.14, 

H  =  (e  A;a-20)A  +  (20-  r«^0)B  -  (^(^+20)  C 
b  b  b 


(4.12) 

(4.13) 


(4.14) 


(4.15) 


(4.16) 
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also,  at  X  =  0 


2  4 

E  =  a'n  T  =  -A/c  -  Bk  +  D 
o  o 


(4.17) 


while  at  X  =  d, 


E,  =  -3a'n^T  ^  +  4ci'n^T  =  -A/ce^^*^  -  B/ce  +  CJd  +  D 
d  o  o  d 


(4.18) 


Then,  equations  4.12,  4.15,  4.17  and  4.16  constitute  the 
boundary  condition^:.  a  set  of  four  simultaneous  linear  al¬ 
gebraic  equations  in  A,  B,  C  and  D; 


I  =  (1-^)  A  +  (1+8)  B  -8C  +  D 

o  <* 


(4.12) 


E  =5 

o 


-  KB 


(4.17) 


_  ad^  -ad 

E„  =  -KB  A  -  /ce  B 

D 


+  dac  +  D  (4,18) 


H  =  'J:  «;a  -  28)A  -  /ca  -  28  )  B  -  (^  a  +  28)C  (4.16) 

b  o  b 


Solving  these  equations  for  A,  B,  C,  and  D,  we  obtain. 


for  I  =  0, 
o 


A  =  j  fiiiliil  [.E  d+B+zce”^"^)  +  E,(1+8+/c)] 
A  /c  o  a 


+  H[da(l+8+K.)  +  ^;8(l-e  )]  } 


(4.19) 


B  =  “  -  E^(1-8+/c)] 

ord 

-  H[da(l-8+/c)  +  «:3(l-e  )])  (4.20) 

C  =  7  Hk{(1-3>/c)  d-e"^^)  -  (1+3+/C)  (1-3^^)  }  (4.21) 

A 

F  =  E  +  T/c  ( — [E  (/c  sinh  cjd-3)  +  3e  ] 
o  a  /c  o  d 

+  H  [23  (da  +  ic  sinh  ad)])  (4.22) 

=  E  +IC  (A+B) 
o 


where 

A  =  (1+k;)  4  3  {3  (1-cosh  ad)  -  (l+/c)sinh  ad)  (4.23) 

The  suitability  of  the  Hamaker  solution  for  this  kind  of 
boundary  conditions  is  one  of  its  strong  features.  In  contrast, 
there  is  no  obvious  way  in  which  the  exact  solution  can  be 
utilized  except  when  the  four  necessary  boundary  conditions  are 
all  given  at  one  point. 
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5.0  SAMPLE  PREPARATION  AND  CHARACTERIZATION 


Samples  of  several  materials  were  prepared  for  emissivity 
measurements.  The  materials  and  cheir  sources  and  average  par¬ 
ticle  size  were; 


Material 

Powder 

Particle  Size 

Fabricator 

Al^Os 

Linde  A 

0.3  micron 

Lexington 

Laboratories 

MgO 

1  micron 

Avco  RAD 

SrTiO^ 

TAM  C.P. 

-325  mesh 

Lexington 

Laboratories 

SiO^ 

Avco  RAD 
&  Corning 

CaF^ 

Baker  CP 

Lexington 

Laboratories 

Initial  attempts  were  made  to  prepare  all  of  the 
materials  except  SiO^  by  isostatic  pressing  and  sintering.  This 
procedure  resulted  in  complete  success  with  the  Al^O^/  partial 
success  with  SrTiO^  and  CaF^*  Other  methods  were  used  to  pre¬ 
pare  the  remaining  samples. 

Measurements  of  density  were  performed  by  both  point 
count  on  a  polished  section  (counting  the  pores  falling  under 
the  intersections  of  a  uniform  grid)  and  by  immersion  in  xylene. 
Xylene  was  chosen  because  of  its  low  surface  tension  relative 
to  the  solids  to  be  immersed.  Weights  of  the  sample  were  taken 
(a)  dry,  (b)  immersed  and  (c)  saturated  but  not  immersed.  The 
measurements  have  been  generally  reproducible  within  1%. 

Pore  radius  was  computed  from  the  measured  average  length 
of  line  intersected  by  the  pores.  The  relation  used  to  obtain 
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the  value  of  the  pore  radius  was 


r 


(5.1) 


where  r  =  average  pore  radius  and  &  ~  average  length  of  line 
intercepted  by  pore. 

5.1  AI2O3 


5.1.1  Preparation 

Five  samples  of  Al^O^  were  prepared  by  isostatic 
pressing  at  8000  psi.  An  addition  of  20  wt  %  distilled  water 
was  made  to  aid  green  strength.  The  samples,  as  pressed,  were 
disks  approximately  5  inches  in  diameter  by  2  inches  in  thick¬ 
ness;  extra  material  was  allowed  in  order  that  the  surface  be 
characteristic  of  the  bulk  material. 

Firing  was  done  in  air  using  a  furnace  with  a 
program  controller.  Variations  of  time  and  temperature  were 
used  to  control  the  density.  Where  necessary  samples  were 
fired  more  than  once  to  adjust  the  density  to  a  desired  value. 


5.1.2 

follows 

Characterization 

The  measured 

{%  of  theoretical) : 

density  of  the  Al^O^ 

samples  is  as 

Sample 

Point  Count 

Immersion 

Figure 

A.2-1 

88.1 

95.8 

5.1 

Ai-2 

77.4 

86.0 

5.2 

Ai-3 

96.4 

95.8 

5.3 

Ai-4 

85,5 

90.0 

5.4 

Ai-6 

77.0 

87.6 

5.5 

Greater 

discrepancies  were  observed  between  the 

less  dense 

samples 

for  two  reasons;  (1)  low  density  and  fine  grain  size 
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5 . 2  MgO 

5.2.1  Preparation 

Controlled  ^t  pressing,  performed  at  Avco  RAD, 
was  done  in  graphite  dies.  A  measured  charge  was  placed  in 
the  die,  then  pressure  and  temperature  increased  until  a  pre¬ 
determined  ram  travel  was  observed. 

Machining  of  the  finished  samples  v'as  done  dry 
using  diamond  tools.  Special  care  was  observed  with  the  lower 
density  sample  because  of  the  tendency  of  fine  particle  MgO 
to  react  with  atmospheric  CO^  and  H^O. 

5.2.2  Characterization 

Density  measurements  provided  the  following 
values  (%  of  theoretical) : 

Sample  Point  Count  immersion  Figure 

Mg-3  88»0  94.0  5.6 

Mg-5  90.0  97.5  5.7 

The  fine  grain  size  of  the  less  dense  material 
generated  exceptional  difficulty  in  producing  a  polished  section 
that  could  be  photographed. 
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was 


The  average  pore  radius  observed  for  each  sample 


Sample 
flg-3 
Mg- 5 


Pore  Radius  (microns) 


2.13 

1.78 


5.3  SrTiO^ 

5.3.1  Preparation 

Isostatic  pressing  and  sintering  was  used  to 
obtain  a  high  density  sample  (Sr-1) ;  slip  casting  was  used  to 
obtain  a  lower  density  sample.  Both  samples  were  fired  in 
air  using  slow  heating  rates. 

Samples  were  machined  wet  using  diamond  tools. 

5.3.2  Characterization 

Sample  densities  were  observed  as  follows  {%  of 

theoretical) : 

Sample  Point  Count  Immersion  Figure 

Sr-1  92.9  98-7  5.8 

Sr-2  (cast)  *  63.  64.  5.9 

Average  pore  radius  for  the  two  samples  was 
Sample  Pore  Radius  (microns) 

Sr-1  2.47 

Sr-2  (cast)  3.00 


5.4  siO^ 


5.4.1  Preparation 

Two  sample  preparation  techniques  were  used;  a 
higher  density  sample  was  hot  pressed  by  Avco,  and  two  of  lower 
densities  were  slip  cast  bv  Corning,  and  fired  to  a  terminal 


*  Because  of  the  difficulty  in  polishing  Sr-2  (cast) ,  the 
value  .63  was  used;  it  was  obtained  both  on  a  xylene  im¬ 
mersion  measurement  and  by  mechanical  measurement  of  the 
weight  and  volume  of  a  regular  cylinder  of  the  material. 


^0 


density.  In  both  cases,  the  maximum  temperature  was  low  enough 
to  cause  little  or  no  devitrification. 

The  hot  pressed  sample,  however,  was  observed 
to  be  over  99%  dense.  As  a  result,  it  was  not  considered  a 
semi-transparent  material  with  pores  as  scattering  centers,  but 
a  nearly  transparent  glass.  The  two  samples  from  Corning  were 
measured  and  the  resulting  data  used. 

Calculated  values  for  emissivity  of  fused  silica 
were  computed  on  the  basis  of  the  microstructure  of  the  Corning 
samples. 

Machining  was  done  wet  with  diamond  tools. 

5.4.2  Characterization 

Observed  values  of  density  were  (%  of  theoretical) 


Sample 

Point  Count 

Immersion 

Figure 

FS-3 

65.9 

84.3 

5.10 

PS-5 

71.7 

87.0 

5.11 

FS-7 

100.0 

99.5 

5.12 

Average  pore  radius  for  the  two  samples  was: 
Sample  Average  Pore  Radius  (microns) 

FS-3  1.28 

FS-5  3.06 

5.5  Cr.F^ 


5.5-1  Preparation 

Efforts  at  preparing  a  satisfactory  sample  were 
thwarted  by  adsorption  of  water  on  the  surface  of  the  powder; 
during  firing  the  water  reacted  to  form  a  complex  hydrate  and 
finally  CaO.  Attempts  at  preparing  the  powder  and  firing  the 
slug  in  an  atmosphere  of  dry,  CO^ -free  nitrogen  were  partially 
successful.  Two  samples,  free  of  CaO,  were  prepared;  but  failed 
in  machining. 
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Fig.  5.1  Ai-1  532X.  Point  Count  Density  =  88.1%  of 
Theoretical  Density. 


Fig.  5.2  A£*2  532X.  Point  Count  Density  =  77.4%  of 

Theoretical  Density. 
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Fig.  5.3  Ai-3  532x.  Point  Count  Density  =  96.4%  of 
Theoretical  Density. 


Fig.  5.4  A^-4  532x.  Point  Count  Density  =  85.5%  of 

Density . 
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4  »» 


Fig.  5.5  Ai-6  532X.  Point  Count  Density  =  77.0%  of 
Theoretical  Density. 


Fig.  5.6  Mg-3  532X.  Point  Count  Density  =  88.0%  of 
Theoretical  Density. 
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Fig.  5.7  ^^-5  532X.  Point  Count  Density  =  90.0%  of 

Theoretical  Density. 


Fig.  13. &  Sr-1  1064X .  Point  Count  Density  =  92.9%  of 
Theoretical  Density. 
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Fig.  5.9  Sr-2  1064X.  Immersion  Density  =  63%  of  Theoretical 
Density. 


Fig.  5.10  FS-3  532X.  Point  Count  Density  =  65.9%  of 
Theoretical  Density. 
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Fig.  5.11  FS~5  532Xo  Point  Count  Density  =  71.7%  of 
Theoretical  Density. 


Fig.  5.12  FS-7  532X.  Polarized  Incident  Light  Shows 

Crystalline  Phase.  Immersion  Density  =  99.5%  of 
Theoretical  Density. 
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6.0  SINGLE  CRYSTAL  ABSORPTIVITY 


Transmission  of  various  single  crystals  was  measured  to 
determine  the  adDsorptivity ,  or,  as  a  fxinction  of  wavelength  and 
temperature.  Measurements  were  made  at  temperatures  up  to 
1200®C  and  at  wavelengths  as  long  as  15  microns  or  as  permitted 
by  material  properties.  In  the  calculation  of  o',  correction 
was  made  for  surface  reflections  by  the  following  formula: 


(1-p) 


2  -at 
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(6.1) 


where 


The  value  of 
if  the  index 
The  relation 


I  =  transmitted  energy 
I^  =  incident  energy 
t  =  thickness  of  sample  (cm) 
e  =  natural  logarithm  base 
p  =  surface  reflectivity 
p  at  a  particular  wavelength  may  be  calculated 
of  refraction,  n^,  is  known  at  that  wavelength, 
is 


P 


2 

(1-n^) 

- —2 

(1+a^) 


(6.2) 


6.1  Experimental  Technique 

A  small  high  temperature  furnace  with  a  rectangular 

cross-section  was  used  to  heat  the  samples.  Special  optics 

are  required  for  use  with  the  spectrophotometer.  Details  have 

(3) 

been  presented  previcusly. 

6 . 2  Samples 

Single  crystals  of  Al^O^#  MgO,  SrTiO^  and  CaF^,  were 
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measured  as  well  as  s^unples  of  Si02  glass.  Isothermal  emis- 

sivity  measurements  were  done  on  most  of  the  same  group  of 

materials.  Sources  of  the  materials  were 

Al^O^  Linde 

MgO  Semi-Elements 

SiO„  A.  D.  Jones  Optical 

£ 

SrTiO  TAM  Division 

National  Lead 

CaF-  Semi-Elements 

6.3  Results 

Figures  6.1  through  6.5  show  absorption  coefficients 
as  a  function  of  temperature  and  wavelength. 
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Fig.  6.1  Absorption  coefficient  of  Single  Crystal  Al^O^* 
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Pig.  6.2  Absorption  coefficient  of  Single  Crystal  MgO. 
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Absorption  coefficient  of  Single  Crystal  SrTiO 
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Fig.  6.5  Absorption  coefficient  of  Si02  Glass- 
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7.0  ISOTHERlilAL  EMISSIVITY  AND  CALIBRATION  MEASUREMENTS 

Measurements  of  isothermal  emissivity  were  made  on  well 
characterized  samples  and  a  check  of  the  effectiveness  of  the 
blackbody  source  was  performed  using  a  special  blackbody  as  a 
reference. 

7.1  Isothermal  Emissivity 

7.1.1  Measurements 

Measurements  of  emissivity  were  performec  using 
the  same  equipment  as  previously  described.^  Basically  it 
consists  of  two  disk  shaped  scunples  adjoining  each  other  within 
a  furnace  with  optical  access  through  a  water  cooled  port. 
Transfer  optics  provide  two  beams  to  a  recording  spectrophotom¬ 
eter,  one  the  reference  beam  from  a  cylindrical  groove  between 
the  sar^nles  that  provides  a  blackbody;  the  other  from  the  cylin¬ 
drical  surface  of  the  sample  as  it  rotates  past  the  port.  Thus 
the  spectrophotometer  records  emissivity  of  the  sample  directly 
with  no  correction  for  temperature. 

7.1.2  Results  of  isothermal  measurements 

Figures  7.1  through  7.29  present  the  results 
of  measurements.  Figures  7.30  through  7.36  present  comparison 
calculations  made  during  the  present  contract  and  measurements 
made  on  the  previous  contract. 

7.2  Blackbody  Calibration 

Emissivity  measurements  made  by  direct  comparison 
of  the  radiation  emitted  from  a  sample  surface  and  a  reference 
blackbody  can  only  be  as  good  as  the  blackbody  reference.  Two 
general  requirements  must  be  met:  (1)  the  blackbody  must  have 
an  emissivity  near  unity  and  (2)  the  blackbody  must  be  at  the 
same  temperature  as  the  seimple. 
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In  the  Lexington  Laboratories  design,  these  require¬ 
ments  are  met  with  a  rotating  wedge  which  is  actually  part  of 
the  sample  and  thus  is  at  the  same  temperature.  The  slot  sur¬ 
face  is  blackened  with  a  silicon  carbide  coating  having  an 
emissivity  near  0.9  over  most  of  the  wavelengths  of  interest 
and  above  0.7  at  all  wavelength,  wi  h  this  surface,  the  geom¬ 
etry  is  such  as  to  give  an  emissivit  above  0.995  over  all 
wavelengths  of  interest. 

7.2.1  Apparatus 

A  blackbody  furnace  has  been  constructed  and 
the  essential  details  of  the  unit  are  shown  in  Fig.  7.37. 

The  unit  is  based  on  a  design  used  at  the  National  Bureau  of 
Standards. By  their  calculations,  the  calculated  emissivity 
was  0.999,  although  the  actual  emissivity  is  probably  some¬ 
what  less. 

The  only  changes  in  the  NBS  design  are:  (a) 
the  use  of  Kanthal  A-1  instead  of  Pt-20%Rh,  (b)  machining  the 
core  instead  of  casting,  and  (c)  use  of  air  cooling  instead 
of  water. 

The  furnace  is  introduced  to  the  optical  system 
of  the  spectrophotometer  by  rotating  one  of  the  mirrors  in  the 
external  entrance  optics  (See  Fig.  7.38).  By  rotating  this 
mirror,  the  system  may  be  adjusted  to  the  proper  energy  source 
without  introducing  any  auxiliary  mirrors  but  still  using  the 
full  complement  of  mirrors.  Introduction  of  auxiliary  mirrors, 
exclusion  of  any  mirrors,  or  changes  in  reflection  angle  would 
tend  to  introduce  differences  in  the  beam  path.  Such  differ¬ 
ences  would  prevent  useful  calibration. 

7.2.2  Calculations 

In  emissivity  measurements  and  comparison  tests 
using  a  separate  blackbody,  it  is  essential  to  have  well-defined 
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temperatures.  The  energy  emitted  by  a  blackbcxJy  as  a  funct  on 
of  temperature  is  given  by  Planck's  rac'iition  law: 


-5  1 


C,/AT 
e  -1 


(7.1 


For  considerirg  the  effect  of  small  temperat  re 
differences,  Wien's  radiation  law  is  a  good  approximation. 


^7\,T  ^  e-C2AT 


Then,  if  the  reference  blackbody  and  the  san  le 
are  at  slightly  different  temperatures. 
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where  =  1.435  when  A  is  given  in  centimeters  and  T  in 
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7.2.3  Results  auid  discussion 

Inspection  of  equation  7.5  shows  that  the  effect 
of  temperature  differences  is  particularly  important  at  the 
shorter  wavelengths.  That  is,  at  1000®C,  a  difference  in  tem¬ 
perature  of  5®  between  the  standard  and  the  reference  leads 
to  a  difference  of  energy  emittec.  by  the  blackbody  of  about 
5%  at  1  micron  and  only  .5%  at  ten  microns.  Thus,  in  comparing 
blackbodies,  measiirements  at  the  snorter  wavelengths  are  of 
particular  significance. 

Results  of  one  set  of  measurements  are  illustrated 
in  Fig.  7.39.  The  blackbody  furnace  was  at  a  temperature  of 
100e®C  and  stabilized;  the  rotating  slot  was  stabilized  at  a 
temperature  of  10'^3®C.  The  ratio  of  measured  energy  emitted 
was  within  less  than  1%  of  the  calculated  value  over  the  entire 
range  of  wavelengths  and  no  systematic  deviation  was  observed, 
except  perhaps  at  the  shorter  wavelengths  where  small  variations 
in  temperature  assume  greater  significance.  We  interpret  these 
results  as  confirming  the  blackness  of  our  rotating  wedge  con¬ 
figuration,  and  illustrating  the  need  for  precise  temperature 
measurement  and  control  when  an  exterior  blackbody  source  is 
used. 

In  the  case  of  these  measurements  we  also  con¬ 
firmed  previous  observations  of  the  need  for  care  in  alignment, 
particularly  at  the  longer  wavelengths  when  the  spectrometer 
slit  is  wide  open.  A  measurement  under  similar  conditions  for 
those  shown  in  Fig.  7.39  is  illustrated  in  Fig.  7.40.  Here 
the  agreement  between  blackbody  furnace  and  rotating  wedge  is 
good  at  the  shorter  wavelengths,  but  decreased  energy  is  ob¬ 
tained  from  the  rotating  wedge  at  wavelengths  above  about  11 
microns  because  of  misalignment. 
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Thus  the  present  results  confirm  our  previous 
experience  that  experimental  measurements  are  reproducible 
within  +  1%.  However,  the  possibility  of  minor  misalignment, 
changes  in  optical  path,  or  other  minor  variations,  leads 
us  to  conservatively  claim  no  more  than  a  maximum  possible 
error  of  up  to  5%  for  any  individual  measurement. 
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Pig.  7.7  Measured  vs.  Calculated  Emissivity  of  A1  O  (Ai-3)  590 
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Fig.  7-9  Measured  vs.  Calculated  Emissivity  of  A1  O  (A*.- 3)  1108 
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Pig.  7.14  Measured  vs.  Calculated  Emissivity  of  A1  (Ai-6)  900 
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Fig.  7.15  Measured  vs.  Calculated  Emissivity  of  Al  O  (Ai-6)  1150 
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Fig.  7.21  Measured  vs.  Calculated  Emissivity  of  MgO  (Mg-3)  908 
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Fig..  7.23  Measaired  vs.  Calculated  ESnissivity  of  MgO  (Mg-5)  595 
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Pig.  7.25  Measured  vs.  Calculated  Emissivity  of  MgO  (Mg-5)  1168 
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Fig.  7.26  Measured  vs.  Calculated  Emissivity  of  SrTio  (Sr-1)  605 
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Fig.  7.28  Measured  vs.  Calculated  Emissivity  of  SrTiO  (Sr-2)  573 
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Fig.  7.31  Measured  vs.  Calculated  Emissivity  of  Commercial  A1„0_ (AD-96) 
910°C.  2  ^ 
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7.35  Measured  vs.  Calculated  Emissivity  of  Commercial  Al.,0  (AV-30) 
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Fig.  7.36  Measured  vs.  Calculated  Emissivity  of  Commercial  A1 
(AP-35)  910°C. 


Fig.  7.37  Blackbody  Furnace 
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Fig.  7.40  Ratio  of  Energy  Emitted  by  Blackbody  Furnace  to  t 
Emitted  by  Rotating  Slot  with  Slight  Misalignment 


8.0  COMPARISON  OF  MEASURED  AND  CALCULATED  EMI  SSI VI TIES 


Directly  measured  values  of  the  normal  emissivity  are 
compared  with  values  calculated  from  our  material  char ac ter iza-- 
tion  in  Figs.  7.1  -7.29. 

8.1  Aluminum  Oxide 

Measured  and  calculated  values  for  specially  prepared 
high  purity  aluminum  oxide  are  shown  in  Fig.  7.1-7.15.  In 
the  2-4  micron  range  where  absorption  is  low  and  the  scatter¬ 
ing  coefficient  high  the  emissivity  is  low;  calculated  and 
measured  values  generally  agrie  to  better  than  +0.05.  In  some 
cases  the  calculated  values  are  greater  than,  and  in  some  cases 
less  than,  the  measured  values.  Inasmuch  as  the  sum  of  the 
precision  of  measurement  and  calculation  is  abouc  jt  0.05,  we 
regard  the  agreement  as  good. 

At  wavelengths  of  1  and  5  microns,  where  the  absorp¬ 
tion  coefficient  is  increasing  rapidly,  the  calculated  values 
are  in  all  cases  slightly  or  significantly  lower  than  the  meas¬ 
ured  values.  This  results  from  the  fact  that  the  purity  of 
the  polycrystalline  samples  is  less  than  the  single  crystals, 
and  thus  the  absorption  coefficient  is  less  for  the  single 
crystals,  giving  a  low  value  for  the  calculated  emissivity. 

In  all  cases  the  general  shape  of  the  calculated 
and  measured  curves  are  similar. 

8.2  Silica  Glass 

Measured  and  calculated  values  for  silica  glass 
samples  are  shown  in  Fig.  7.16  -  7.19.  As  with  the  alumina 
samples,  calculated  and  measured  values  are  in  close  agree¬ 
ment  in  the  2-4  micron  range.  Calculated  values  at  1  and  5 
microns  are  lower  than  measured  values,  and  the  general  shape 
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of  t’ e  einispivity  wa,/elength  curve  is  the  same  for  the  calcu¬ 
lated  and  measured  values. 

3.3  Magnesia 

Measured  and  calculated  values  for  hot-pressed  magne- 
sixim  oxide  samples  are  shown  in  Figs.  7.20  -  7.25.  Here  the 
influence  of  the  absorption  coefficient  is  particularly  pro¬ 
nounced.  The  hot-prossed  samples  were  observed  to  darken 
dxiring  testing,  preau?iibly  because  of  carbon  pick-up  in  the 
hot  pressing  operation  (dene  in  graphite  dies) .  As  a  result, 
the  measiired  emissivity  is  several  times  as  large  as  that  cal¬ 
culated;  that  is,  the  single  crystal  absorption  data  are  com¬ 
plete  Ly  unsuitable  for  the  hot-pressed  pclycrystailine  samples. 

8.4  Strontium  Titanate 

Measured  and  calculated  values  for  sintered  strontiiini 
titanate  are  shown  in  Figs,  7.26  -  7.29.  Sample  SR-1  was  sin¬ 
tered  for  a  long  time  at  high  temperatures  and  was  visually 
observed  to  darken.  a  result,  the  measured  emicsivities 
are  sxibst antially  greater  than  those  calculated  from  single 
crystal  data.  Sample  SP-2  (having  higher  porosity)  was  sin¬ 
tered  for  a  short  time,  darkened  less,  and  the  measured  and 
calculated  values  are  in  much  better  agreement. 

Because  of  the  tendency  for  the  absorption  coeffi¬ 
cient  of  strontium  titanate  to  cliange  as  a  function  of  heat 
treatment,  the  emissivity  will  also  change  as  a  function  of 
thermal  history. 

8.5  Commercial  Alumina 

A  comparison  betv/een  the  measured  and  calculated 
emissivities  for  several  samples  of  commercially  prepared 
alumina  are  shown  in  Figs.  7.30  -  7.36.  For  the  pxirest  of 
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these,  A0995,  agreement  between  calculated  and  measured  values 
is  reasonably  good.  For  the  least  pure,  AD-85,  measured  and 
calculated  values  differ  by  a  factor  of  four.  This  results 
primarily  from  the  fact  that  calculated  values  are  based  on 
the  absorption  coefficient  of  a  single  crystal,  while  actual 
absorption  coefficients  of  the  impure  polycrystalline  samples 
increase  with  the  impurity  content. 
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9.0  APPARATUS  AND  METHOD  FOR  MEASURING  NON> ISOTHERMAL 

EMISSIVITY 

The  pxirpose  of  the  apparatus  herein  described  is  to  pro¬ 
duce  a  controlled  temperature  gradient  across  a  disc- shaped 
ceramic  scunple  whose  emissive  properties  are  to  be  measured. 

The  gradient  is  produced  only  in  the  radial  direction,  there¬ 
by  reducing  the  analysis  to  that  of  a  one  dimensional  heat 
transfer  problem.  A  rotating  sample  arrcuigement  is  used  to 
maintain  symmetric  temperature  conditions. 

The  apparatus  operates  inside  the  furnace  which  has  been 
used  for  isothermal  emissivity  mear.urement;  the  latter  has 
been  described  previously.  The  furnace  heating  elements  pro¬ 
duce  either  a  uniform  temperature  thrcwghout  the  sample,  or 
adjust  the  surface  temperature.  An  internal  resistance  ele¬ 
ment  is  used  to  produce  a  positive  gradient  across  the  S"'mple 
while  a  hollow  shaft  i a  used  to  cool  the  samples  with  forced 
air.  The  complete  apparatus  is  shown  in  Figure  9.1. 

The  samples  and  the  internal  heating  element  are  supported 
on  a  one-half  inch  alumina  puaft.  The  alumina  shaft  is  sup¬ 
ported  and  rotated  by  Inconel  shafts  with  conical  points  to 
mate  with  conical  holes  at  each  end  of  the  alximina  shaft. 

The  two  Inconel  shafts  are  rotated  by  an  elec: trie  motor  and 
synchronized  by  a  pair  of  selsyns.  This  arrangement  assures 
minimum  stress  and  vibration  on  the  shaft  and  sample  assembly 
while  giving  a  uniform  rata  of  rotation. 

Provision  has  been  made  for  cooling  the  interior  of  the 
sample  by  air.  A  hollow  alumina  shaft  is  used  in  this  arrange¬ 
ment  wi  h  compressed  air  forced  through  this  shaft.  This 
arrangement  has  not  been  tested. 
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shell  (water  cooled ) 
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Pig.  9.1  Non-Iso thermal  Emissivity  Apparatus 


9.1  Production  of  Radial  Thennal  Gradient 

In  order  to  minimize  the  effects  of  thermal  varia¬ 
tions,  the  apparatus  was  designed  to  produce  a  radial  tempera¬ 
ture  gradient  in  the  sample.  The  gradient  is  obtained  in  a 
sample  which  has  a  width  of  one-half  inch  and  which  is  placed 
between  ceramic  discs  of  the  same  dimensions  and  made  of  a 
material  with  similar  thermal  properties.  The  discs  are  held 
together  by  an  axial  spring  load;  the  resulting  thin  air  gap 
between  adjacent  discs  reduces  eocial  heat  transfer  and  result¬ 
ant  temperature  gradients  along  the  samples  and  support  plates. 

The  interior  heating  element  is  wound  for  a  length 
of  nine  inches  under  the  sample  and  support  plates,  with  the 
samples  located  at  the  center  of  the  winding.  This  arrange¬ 
ment  gives  a  uniform  axial  heat  zone  in  the  region  of  the 
sample  and  support  plate  assembly;  no  axial  gradient  is  pos¬ 
sible  since  there  is  no  variation  in  thermal  properties  in 
the  axial  direction. 

9.2  Internal  Heating  Element 

The  internal  heating  element,  producing  the  non- 
isothermal  condition  consists  of  a  Pt-40%  Rh  alloy  resistance 
element  rotating  with  the  samples.  This  wire  element  is  spi¬ 
rally  wound  around  a  one-inch  O.D.  alumina  core.  The  core 
is  grooved  for  a  ten  turns  per  inch  winding  and  flame  sprayed 
with  alumina  to  cover  the  resistance  wire.  The  outer  surface 
is  ground  true  to  facilitate  assembly  of  the  sample  and  sup¬ 
port  plates  on  the  heating  element.  The  power  for  this  heater 
is  supplied  to  the  rotating  element  through  two  slip-ring 
assemblies  located  on  the  drive  shafts  at  either  end  of  the 
furnace . 
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9.3  Measurement  of  the  Temperature  Gradient 

The  radial  temperatiire  and  temperature  gradient  in 

the  disc  shaped  scunples  are  measured  by  eight  thermocouples 
imbedded  in  the  samples  and  supported  by  the  backup  plates. 
These  thermocouples  are  arranged  in  two  groups  of  four  to 
give  two  sets  of  measurements  for  locations  180°  apart  on 
the  scunple.  Each  set  of  couples  is  positioned  on  a  radial 
spiral  covering  90°  of  sample.  This  minimizes  any  effect 
the  couples  have  on  the  radial  heat  transfer  or  temperature 
distribution  in  the  scunples.  Normally  two  samples  will  be 
run  adjacently  with  one  set  of  thermocouples  in  each. 

The  thermocouple  elements  are  Chrome 1- A lumel  en¬ 
cased  in  an  Inconel  sheath.  The  thermocouples  may  be  wired 
to  indicate  both  absolute  and  relative  temperatures.  Signals 
from  the  thermocouples  are  received  at  one  end  of  the  furnace 
by  a  low  noise  slip  ring  assembly  mounted  on  the  exterior 
drive  shaft.  Measurement  is  done  with  a  Leeds  &  Northrup 
K-3  potentiometer. 

9.4  Operation 

The  seane  optics  used  for  isothermal  emissivity  meas¬ 
urements  are  used  in  this  experiment  with  one  change.  As 
shown  in  Pig.  7.38,  one  mirror  is  rotated  to  receive  the  radia¬ 
tion  from  a  blackbody  furnace.  The  latter  furnace  is  main¬ 
tained  at  a  temperature  near  the  surface  temperature  of  the 
sample . 

9.5  Non-Isothermal  Measurements 

Fig.  9.2  is  a  plot  of  the  ratio  of  the  emissivity  of 
an  AI2O2  body  under  isothermal  conditions  (radiating  surface 
at  895°C)  to  the  emissivity  of  the  same  body  (radiating  sur¬ 
face  at  898° C)  with  a  positive  temperature  gradient  across 
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WAVELENGTH  (MICRONS) 

Fig.  9.2  Comparison  of  Isothermal  and  Non-Isothermal 

Emissivity  of  an  AI2O3  Body  (western  Gold  and 
Platinum  Company,  Al-3-96%,  Si02'-4%)  .  Non- 
Isothermal  Radiating  Surface  at  &98‘’c  -  Thermal 
Gradient  60°C/cm.  Isothermal  Radiating  Surface 
at  895®C. 
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the  sample  of  114°C  (60°C/cm) .  In  both  cases  measurement  of 
emissivity  was  made  relative  to  an  external  blackbody  at  approx 
imately  the  same  temperature  as  the  radiating  surface.  Tempera 
ture  differences  between  blackbodies  cmd  radiating  surfaces 
were  corrected  using  Wien's  law  (Eg.  7.5). 
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10.0  DISCUSSION 


Over  the  wavelength  from  0.5-15  microns  different  processes 
affect  the  emissivity  of  oxide  materials  and  as  a  result  no 
single  analysis  is  suitable  for  the  entire  range  of  material 
behavior.  At  wavelengths  shorter  than  about  6  microns  the 
absorption  coefficient  is  sufficiently  low  that  emissivity  is 
a  volume  process  and  is  strongly  affected  by  the  absorption  co¬ 
efficient  and  scattering  coefficient  of  the  oxide.  At  wave¬ 
lengths  above  about  6  microns  the  oxides  are  essentially  opaque; 
as  a  result,  emissivity  is  a  surface  characteristic  which  is 
affected  by  the  surface  reflectivity  (depending  on  the  index 
of  refraction)  and  surface  roughness,  and  in  general  is  simi¬ 
lar  to  metals  and  other  opaque  solids  as  far  as  the  material 
characteristics  which  need  to  be  defined.  Because  the  reflec¬ 
tivity  is  low,  the  emissivity  is  high  in  the  range  up  to  about 
11  microns.  At  higher  wavelengths,  as  the  reststrahl  frequency 
is  approached,  the  absorption  coefficient  remains  high,  but 
the  reflectivity  increases  and  as  a  result,  the  emissivity 
drops  off. 

Under  conditions  where  oxides  are  opaque,  that  is  at  wave¬ 
lengths  greater  than  6  microns  or  so,  the  process  is  not  sig¬ 
nificantly  different  than  for  other  opaque  solids  and  does 
not  offer  significantl’^  new  problems.  At  wavelengths  shorter 
than  about  6  microns,  the  process  is  quite  dissimilar  from 
that  in  most  materials  which  have  been  studied  previously, 
and  as  a  result,  it  is  in  this  wavelength  range  that  we  have 
concentrated  our  efforts. 

In  the  region  where  a  high  absorption  coefficient  does  not 
overwhelm  the  voliime  process,  the  emissivity  depends  on  the  re¬ 
flectivity  of  the  surface  (a  function  of  the  index  of  refraction) 
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and  on  the  material  constant  $  =  -J  a/(a+2s) .  Thus,  in  order 

o 

to  estimate  the  emissivity,  it  is  necesscury  to  calculate  the 
values  of  the  absorption  coefficient  and  scattering  coeffi¬ 
cient.  Values  of  the  absorption  coefficient  have  been  deter¬ 
mined  from  transmissivity  measurements  with  single  crystals. 
Values  for  the  scattering  coefficient  have  been  calculated  on 
the  basis  of  pore  size,  relative  index  of  refraction,  wave¬ 
length  and  pore  concentration.  The  precision  of  this  calcu¬ 
lation  (which  requires  knowledge  of  the  back  scattering  co¬ 
efficient  as  well  as  total  scattering)  is  not  high,  but  since 
it  enters  as  a  square  root  relationship  in  the  determination 
of  emissivity,  the  precision  of  the  final  emissivity  calcula¬ 
tion  is  better  than  that  for  the  estimation  of  the  scattering 
coefficient. 

A  comparison  of  the  calculated  and  measured  values  of 
emissivity  at  2-4  microns  shows  good  agreement  for  specially 
prepared  high- purity  aluminum  oxide  (within  +0.05  emissivity 
units) .  This  result  is  regarded  as  satisfactory  considering 
that  it  is  determined  entirely  from  material  characteristics 
such  as  microstructure  and  pore  concentration  and  is  inde¬ 
pendent  of  actual  optical  property  measurements.  It  indicates 
the  strong  influence  of  porosity  on  fixing  the  emissivity  - 
the  emissivity  decreases  linearly  with  increasing  porosity. 

It  also,  by  comparison  with  other  saunples  of  alumina,  par- 
ticularly  impure  commercial  materials,  shows  the  great  in¬ 
fluence  of  small  changes  in  composition  and  microstructure 
on  the  resulting  emissivity.  At  a  wavelength  of  3  microns 
impure  commercial  aluminas  have  emissivity  values  as  large 
as  0.5,  while  the  high-purity  and  small  pore  size  samples 
specially  prepared  for  this  study  have  values  as  low  as  0.03. 
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That  is,  for  a  material  such  as  "alumina"  an  order  of  magni¬ 
tude  difference  in  the  emissivity  can  result  from  control 
of  microstructure  and  composition. 

A  comparison  of  the  calculated  and  measured  values  for 
MgO  cuid  strontium  titanate  indicates  the  strong  influence  of 
the  cd>sorption  coefficient  on  experimental  results.  The  MgO 
samples  tested  were  deurkened,  presumcd>ly  resulting  from  carbon 
pick-up  during  hot  pressing  in  graphite  dies,  and  as  a  result 
had  values  of  emissivity  at  3  microns  as  high  as  0.8  compared 
with  calculated  values  of  the  order  of  0.1,  which  would  apply 
to  the  pure  material  witnout  contamination.  Thus,  in  the 
preparation  of  low  emissivity  materials,  contamination  is  of 
significant  importance.  For  strontium  titanate  heat  treatment 
affects  the  absorption  coefficient  (presumably  as  a  result  of 
changed  oxidation  states) .  Thus,  measurements  of  the  emis¬ 
sivity  of  strontium  titanate  and  similar  materials  are  subject 
to  chat^  with  changing  environment,  and  one  must  specify  care¬ 
fully  the  conditions  of  preparation  and  testing  in  order  to 
report  useful  material  data.  This  characteristic  further 
limits  the  utility  of  ted^ulated  emissivity  values  for  oxide 
materials  independent  of  detailed  sample  description. 

In  general  then,  in  the  important  1-5  micron  rcuige  cor¬ 
responding  to  the  range  of  maximum  energy  emission  over  tem¬ 
peratures  from  1000-2500®C  the  emissivity  of  oxide  materials 
can  vary,  or  be  varied,  over  a  wide  range.  For  low  emissivity, 
high  purity,  high  porosity,  and  small  pore  size  is  preferred; 
under  such  conditions  the  emissivity  can  be  well  below  0.1. 

For  high  emissivity,  the  main  requirement  is  to  increase  the 
adasorption  which  can  be  done  by  the  introduction  of  impurities 
or  contauninants  which  give  rise  to  adssorption  in  this  wave¬ 
length  range.  For  materials  such  as  this,  emissivities  of 
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tha  order  of  0.8-0. 9  can  be  achieved.  This  wide  range  in 
values  found  for  oxides  and  a  quantitative  evaluation  of  the 
means  necessary  for  control  of  the  emissivity  is  an  important 
consideration  in  any  application  of  these  materials  at  high 
temperatures  involving  heat  transfer  considerations. 

In  the  measurement  of  isothermal  emissivity  we  have 
compared  our  blackbody  slot  system  (in  which  the  sample  it¬ 
self  corresponds  to  the  reference  standard)  with  an  external 
blackbody  designed  to  have  an  emissivity  greater  than  0.995. 
Experimentally  the  results  with  an  external  blackbody  and  one 
machined  into  the  sample  were  indistinguishable.  The  greater 
ease  in  maintaining  temperature  tiniformity  when  the  blackbody 
slot  is  actually  part  of  a  sample  makes  it  a  preferable  system 
for  experimental  observation. 

In  the  calculation  of  emissivity,  it  is  necessary  to 
know  the  back  scattering  coefficient  for  diffuse  radiation 
and  also  the  absorption  coefficient  for  diffuse  radiation. 

As  a  result,  we  have  analyzed  the  absorption  coefficient  for 
diffuse  radiation  relative  to  the  value  measured  for  plane 
parallel  radiation,  and  also  evaluated  the  scattering  coeffi¬ 
cient  -  particularly  for  systems  in  which  pores  are  the  main 
scattering  centers  and  the  index  of  refraction  of  the  pore 
is  smaller  than  that  of  the  solid  matrix  material.  Under 
these  conditions  the  back  scattering  coefficient  can  be 
evaluated  with  sufficient  precision  to  make  a  useful  estimate 
of  the  resulting  emissivity. 

One  aspect  of  the  problem  has  been  that  we  calculate 
hemispherical  emissivity  while  our  measurements  have  been 
for  normal  emissivity.  After  an  analysis  of  the  transfer 
of  radiation  through  the  boundary  of  the  sample,  we  have 
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evaluated  the  ratio  of  normal  and  hemispherical  emissivity 
and  find  that  it  depends  on  the  relative  index  of  refraction 
but  is  not  far  from  unity  for  most  of  the  oxide  materials  of 
interest. 

Application  of  this  data  to  non- iso thermal  systems  in¬ 
volves  use  of  the  Hamaker  approximation  which  essentially 
linearizes  the  temperature  gradient  in  order  to  obtain  a 
group  of  differential  equations  for  which  an  analytical  solu¬ 
tion  is  available.  Using  numerical  techniques,  we  have  com¬ 
pared  the  Hamaker  and  exact  equations  for  a  series  of  condi¬ 
tions  for  temperature  up  to  3000°C,  for  scattering  coeffi¬ 
cients  up  to  100  cm  and  for  temperature  gradients  up  to 
1000°C/cm.  For  these  conditions  we  find  that  the  Hamaker 
equation  is  in  satisfactory  agreement  with  the  exact  equation 
except  when  values  of  the  lat  ice  conduction  are  much  greater 
than  for  radiant  energy  transfer.  For  ail  practical  applica¬ 
tions  of  oxide  materials  at  temperatures  above  1000°C,  the 
Hamaker  equation  gives  a  satisfactory  approximation  to  the 
exact  equation;  in  fact,  for  values  of  the  absorption  and 
scattering  coefficient  greater  than  unity  it  is  within  1/10 
of  one  per  cent  of  the  exact  equation  for  all  the  conditions 
that  we  have  evaluated. 

One  advantage  of  the  Hamaker  equation  over  the  exact 
solution  is  that  it  is  not  necessary  to  specify  all  of  the 
boundary  conditions  at  one  surface,  and  we  have  extended  the 
Hamaker  equation  to  some  additional  boundary  condit  ions  which 
are  of  greater  utility  than  specifying  all  the  boundary  condi¬ 
tions  at  one  point  in  the  system. 

The  application  of  the  Hamaker  equation  of  most  interest 
is  to  conditions  of  non- isothermal  emissivity.  In  this  case, 
the  effective  emissivity  depends  on  temperature  gradients 
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within  the  body  as  well  as  on  the  surface  temperature  itself. 

We  hav^j  designed  and  tested  an  apparatus  which  allows  us  to 
measure  the  energy  emitted  from  a  surface  of  a  sample  con* 
tailing  a  known  temperature  gradient  in  comparison  with  that 
of  a  blackbody  steuideurd.  Test  measurements  indicate  the 
energy  emitted  from  aliiminum  oxide  cereunics  does  strongly 
depend  on  the  temperature  gradient  as  well  as  on  the  surface 
temperature.  With  a  temperature  gradient  of  BO^C/cm  the 
effective  emissivity  for  the  same  surface  temperature  is 
experimentally  found  to  increase  about  22%  at  1000°C.  In 
subsequent  work  we  pleui  to  carry  out  additional  experimental 
measurements  and  compcure  these  directly  with  calculations 
based  on  sample  chcuracteristics.  As  for  the  case  of  isothermal 
emissivity  it  is  essential  to  have  a  detailed  sample  charac* 
terization  to  make  useful  analyses  cmd  predictions. 
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APPENDIX  I 


VALUES  FOR  INDICES  OF  REFRACTION  USED  IN  EMISSIVITY 

CALCULATIONS* 


X  (u) 

MATERIAL 

0.5/<. 

1 

2 

3 

4 

5 

8  11 

Ai^Oa 

1,78 

1.75 

1.74 

1.72 

1.65 

1.60 

(1.33) 

SrTiO^ 

2,48 

2.31 

2.27 

2.23 

2.16 

2.11 

(1.85) 

9606 

1,45 

1.45 

1.43 

1.41 

(1.35) 

(1.30) 

9608 

1,45 

1.45 

1.43 

1.41 

(1.35) 

(1.30) 

MgO 

1.75 

1.72 

1.70 

1.68 

1.65 

1.55 

(1.45) 

CaF2 

1.43 

1.42 

1.41 

1.41 

1.40 

1.39 

1.35  1.28 

Si02 

1.45 

1.45 

1.43 

1.41 

(1.35) 

(1.30) 

Parentheses  indicate  extrapolated  values. 


*  from  Wolfe  and  Ballard ^ Proc .  Inst.  Radio  Engineers  £7,  1540 
(1959) 
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APPENDIX  II 


SAMPLE  CALCULATION  OF  EMISSIVITY 


II .1  Calculation  of  a  and  a  for  Ai-6  (3u)  at  904°C 

The  calculation  of  emissivity  for  Ai-6  (3|i)  at  904°C  is 
typical  of  the  method  used  for  obtaining  all  values  of  calcu 
la ted  emissivity.  Prom  the  microstructure  of  Ai-6  in  Fig. 
5.5  values  of  P  (pore  volume  fraction)  and  r  (average  pore 
radius)  are  obtained. 

The  absorption  coefficient  for  single  crystal  ^120^  was 
calculated  from  the  ratio  of  energy  transmitted  through  an 
optically  polished  section  0.313  cm  thick  to  that  incident 
upon  it  and  corrected  for  the  surface  reflectance  at  the 
polished  surface  from  the  following  relationship: 


(1  “  p)^ 


-Q-t 

e 


(II. 1) 


where : 


I  * 

I  = 
o 

a  m 
t  * 
The 
index  of 


transmitted  energy 
incident  energy 
absorption  coefficient 
thickness  of  section  in  cm. 

value  of  surface  reflection  is  calculated  from  the 
refraction: 


^  (1  +  n^) ^ 

A 


(II. 2) 


The  value  of  the  absorption  coefficient  a  may  then  be 
determined  by  solving  Eg.  (II. 1) 

a  .  -  iS_Lio - ~  P*  ^  ]  (II. 3) 
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For  Al-0,  at  (3  \x)  and  904®C 

4b  W 


o 


- ^  *  1.155  (at  Tv  =  3  n  =  1.72) 

(1-p) 

t  =  .313  cm. 


When  II. 3  is  evaluated  for  the  stated  conditions 

o  =  .215  cm  ^ 

Equation  2.16  states: 

a  *  2  O'  (2.16) 

Where  a  is  the  absorption  coefficient  for  a  polycrystalline 
material  with  pores  as  scattering  centers,  therefore: 

a  *  .430  cm"^ 

II. 2  Calculation  of  s  for  A£-6  (3n)  at  904° C 

The  scattering  factor  s  is  defined  by  the  expression: 

s  -= — 7 —  K  — - —  (2.44) 

4  r 

P  *  pore  volxoKie  fraction 

r  =  average  pore  radius 

K  «  scattering  factor  from  Fig.  2.3 

4nr  Im-ll _ 

^  *  Tvm  (II. 4) 

where : 

r  »  average  pore  radius  (1.68p) 

1 
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For  Ai-6  (3m)  at  904°C,  P  -  0.23  and  K  -  1.70 
From  2.44  then 

s  ■  1.751  X  lO^can  ^ 

II. 3  Evaluation  of  0^  from  values  of  a  and  s  for  Ai-6  (3m) 
at  904  C 


8 


C 


I  ^ 

-»•  V  a  +  2s 


(II. 5) 


where : 


a  is  the  absorption  coefficient  for  poly¬ 
crystalline  materials  with  pores  as 
scattering  centers.  (.430  cm  ) 

s  is  the  scattering  coefficient  (1.751  x  10 

-1, 

cm  ) 

0  is  a  material  constant 


o  .  /4.30  X  10  cm _ 

^  *  ■*■  /  3  _i 

V  4.30  X  10  cm  +  3.502  X  10  cm 


0. 


..rr. 


23  X  10 


-4 


0Q  *  1.11  X  lO"^ 

II. 4  Calculations  of  the  reflectivity,  p  and  the  emissivity, 
e  from  values  of  0^,  and  p^ 

Values  of  p^  and  p^  are  found  in  Table  2.1 
for  a  range  of  refractive  indices.  Values  between  were  inter¬ 
polated  . 

The  value  of  e  is  derived  from  equation  2.1. 
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Substituting  in  2.1,  we  find: 

p  =  0.938  and 
(1  -  p)  =  e 
€  =  .062 
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(1  +  p^)  +  23q  (pQ  +  p^) (1  -p^) 
+  +  23q(1  -p^^) 


.  constant 

itivity  of  diffuse  energy 
he  sample 

:tivity  of  diffuse  energy 
i  the  sample. 


numerically: 

(1  -p.^)  =  .506 
p  =  .938 

23q  =  2.22  X  10“ 
3q^  *  1.23  X  10“ 

(pO  +  Pi)=  -825 

(1  -  p.  -  2Pq)  *  .053 


